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Abstract 


Ybarra, Gary A. High Resolution Target Range Estimation in Inhomogeneous 
Media Using Millimeterwave Radar*. (Under the direction of Sasan H. Ardalan.) 

Two narrow-band radar systems are developed for high resolution target 
range estimation in inhomogeneous media. They are reformulations of two 
presently existing systems such that high resolution target range estimates may 
be achieved despite the use of narrow bandwidth radar pulses. A double 
sideband suppressed carrier radar technique originally derived in 1962, and 
later abandoned due to its inability to accurately measure target range in the 
presence of an interfering reflection, is rederived to incorporate the presence of 
an interfering reflection. The new derivation shows that the interfering reflection 
causes a periodic perturbation in the measured phase response. A high 
resolution spectral estimation technique is used to extract the period of this 
perturbation leading to accurate target range estimates independent of the 
signal-to-interference ratio. A non-linear optimal signal processing algorithm is 
derived for a frequency-stepped continuous wave radar system. The resolution 
enhancement offered by optimal signal processing of the data over the 
conventional Fourier Transform technique is clearly demonstrated using 
measured radar data. A method for modelling plane wave propagation in 
inhomogeneous media based on transmission line theory is derived and 
studied. Several simulation results including measurement of non-uniform 
electron plasma densities that develop near the heat tiles of a space re-entry 
vehicle are presented which verify the validity of the model. 


* This work was supported in part by the National Aeronautics and Space Administration Langley 
Research Center under Grant NAG-1-1219. 
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Chapter 1 
Introduction 


This dissertation presents the reformulation of two presently existing radar 
systems such that narrow-band pulses may be utilized while achieving high 
resolution target ranging. Conventional pulsed radar achieves high resolution 
by transmitting large bandwidth pulses. If the medium in which the target is 
embedded has a frequency dependent inhomogeneity, then a wide-band pulse 
would be smeared in time due to dispersion because different frequency 
components of the pulse will travel at different velocities. The detection process 
used to measure the target range from the received pulses will produce 
significant error. Therefore, high resolution is not achievable using conventional 
pulsed radar if the propagation path contains any medium with a frequency 
dependent inhomogeneity. 

The two radar systems developed in this dissertation are a frequency- 
stepped double sideband suppressed carrier (DSBSC) system and a frequency- 
stepped continuous wave (CW) system. Specifically, the stepped DSBSC 
process is derived to incorporate the presence of an interfering reflection. The 
derivation shows that the interfering reflection causes a periodic perturbation in 
the measured phase sequence. It is demonstrated that by applying a high 
resolution spectral estimation technique to these measurements, accurate target 
range estimates may be detected. For the stepped CW technique, a new optimal 
signal processing algorithm is derived for extracting highly resolved target range 
estimates. The range resolution achievable using the new optimal processing 
algorithm is clearly demonstrated using measured radar data. 


In order to simulate the performance of the two radar systems, a model for 
electromagnetic (EM) plane wave propagation in inhomogeneous media is 
derived. Simulation results using this model are presented in the context of 
measuring non-uniform electron plasma densities that develop near the thermal 
protection heat tiles of a space re-entry vehicle. However, the model is much 
more general and can be used to simulate accurately EM propagation in any 
inhomogeneous medium. 

1.1 Microwave and Millimeterwave Radar 


There presently exist several schemes which utilize electromagnetic (EM) 
energy to determine the distance from an antenna to a reflective target. The 
process of detecting a target and its range using EM energy is referred to as 
radar, an acronym for Radio Detection and Ranging. The primary objectives of 
a radar system are to determine the presence of a target and also its range, or 
distance from some reference point to the target. Other target properties that are 
capable of being detected by certain radar systems are its location, velocity, and 
physical dimensions [Collins, Skolnik]. The radar systems developed in this 
dissertation are specifically designed for estimating target range. However, with 
further signal signal processing, target velocity can be extracted from the data 
produced by either system. This extension remains as future research and is not 
addressed in this dissertation. 

There are two main performance measures of a target range estimate: 
accuracy and resolution [Artech84, Wehner, Bird]. Accuracy refers to the 
difference between the true and estimated values of the range. Resolution is the 
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ability to detect two or more distinct scatterers when they are at approximately 
the same range. These two performance measures of range estimation are 
often confused and careful attention should be given to their difference. 

The electromagnetic spectrum contains a very wide range of frequencies 
that are currently used for radar measurements. Radars can be designed and 
developed at any frequency at which EM energy can be generated and 
controlled. Most radar systems utilize frequencies in the microwave and 
millimeterwave regions. The frequency ranges associated with microwaves and 
millimeterwaves overlap. Together they span the EM spectrum from 
approximately 300 MHz to 300 GHz. The radars presented in this dissertation 
are designed to operate within this frequency range. 

A radar system is typically designed to achieve superior performance for a 
specific set of measurement circumstances. For example, many applications 
require a long detection range [Artech 84]. In order to detect targets at long 
ranges, careful selection of the radar operating frequency is critical to avoid 
atmospheric absorption. Examples of other measurement circumstances which 
may dictate constraints on system design are need for covertness, tracking 
ability, classification or identification, adverse weather performance, and volume 
search [IEEE]. The radars developed in this dissertation do not consider these 
constraints. They are developed solely for the general purpose of producing 
accurate and highly resolved target range estimates. Specific system details 
that may be required for a particular application remain the responsibility of the 
specific system designer. 
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Traditional pulsed (constant frequency) radars rely on a narrow 
pulsewidth to provide the range resolution Ar, which is directly proportional to 

the pulsewidth [Wehner] 

Ar * ctp (1 .1 .1 ) 

where c is the velocity of EM propagation and tp is the time duration of the radar 

pulse. A popular high resolution variation of traditional pulsed radar is 
frequency-chirp pulsed radar which has a range resolution Ar s 

Ar d -2-1) 

where p is the bandwidth of the frequency-chirped pulse. Frequency chirping is 
achieved by frequency modulating the transmitter pulse. The range resolution 
objective of the radar systems developed in this dissertation is on the order of 5 
mm. In order to achieve this resolution using either traditional pulsed radar or its 
high resolution chirp variation would require a pulse with a minimum bandwidth 
of 30 GHz. if the target is embedded in a frequency dependent inhomogeneity, 
such a pulse would be smeared in time due to the dispersive medium and the 
range estimate will be highly erroneous. Hence, an alternative approach to the 
pulsed radar technique is required. 
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1 .2 The Microwave Reflectometer Ionization Sensor 


The primary motivation for this study has been a prospective NASA 
shuttle mission called the Aeroassist Flight Experiment (AFE). In this 
experiment, a small spacecraft is to be deployed by the shuttle and accelerated 
to atmospheric entry velocity. The unmanned craft will then be flown through the 
earth's atmosphere. During the 600s aeropass, the high temperature in the 
proximity of the non-ablating heat tiles will generate a dynamic electron plasma. 
The aerodynamic friction created will slow the craft, providing an aerobraking 
action. Following the aeropass, the AFE vehicle will attain a low earth orbit to 
be retrieved by the shuttle. One purpose of this experiment is to gather 
information for the development of aerobraking technology, which requires an 
accurate measurement of the electron plasma density profile. The plasma 
profile data will be used to confirm or improve the flow field predictions made by 
computational fluid dynamics. The system which will perform the plasma density 
measurements is called the Microwave Reflectometer Ionization Sensor (MRIS). 
The MRIS antennas will be mounted behind the stagnation region tile where the 
plasma flow field is expected to be laminar. Polarization is not an issue because 
the plasma density is expected to be constant in the transverse plane. 

There exist several predictions of the plasma density profiles that are to 
be measured. One such prediction, based on Computational Fluid Dynamics 
(CFD), is considered to be one of the better predictions presently available 
[Ybarra2]. An example of a CFD predicted profile is shown in Fig. 1 .2.1 . 

The distance values shown are referenced to the outer surface of the heat 
tiles of the AFE vehicle. Distances referenced from the heat tiles are referred to 
as standoff distances. 
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Electron plasma may be characterized by its relative dielectric constant 
[Chen] given by 


e r = 1 - 


q 2 N 

meow 2 (1-P4 

(0 


(1.2.3) 


where q is the electron charge, N is the electron density, m is electron mass, £q 
is the permittivity of free space, to is the angular frequency of propagating 
electromagnetic (EM) energy, and v is the plasma collision frequency. 

CFD ELECTRON DENSITY PROFILE PREDICTION (87 km altitude) 

3 

N (e' /cm ) 



Figure 1.2.1 Computational fluid dynamics (CFD) prediction of the electron 
density profile at 87 km altitude. 
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(1-2.4) 



where the EM frequency f is in kHz, and the electron density N is in e’/cm 3 . 
Unlike other materials, the group velocity v g of EM energy in a plasma is the free* 

space velocity c scaled by yje T . 


Vg = C {Ex = C 



(1 .2.5) 


There exists a critical frequency f„ for a given electron density which occurs 
when e r = 0. At this point in the plasma, called the turning point , the plasma 

becomes a perfect reflector. That is, at the turning point, the reflection coefficient 
becomes T = -1, and interestingly the group velocity goes to zero. For 
frequencies much greater than f CT , £ r ®1 and the plasma behaves similar to free 
space. For Frequencies less than f CT , £ r < 0 the plasma behaves similar to a 

waveguide below its cut-off frequency. That is, only evanescent modes are 
present and no propagation occurs. As an example of a plasma density 

distance measurement, consider the CFD profile in figure 1 .2.1 . Suppose an EM 
wave is launched at a frequency of 140 GHz. The critical density, N c , at this 

frequency may be calculated from (1.2.4) to be 2.42*10 14 e'/cm 3 and is located 
at a standoff distance of 6 cm. The EM wave will penetrate into the plasma until 
it is reflected at the turning point and returns as an echo. If the round-trip travel 
time is measured and an average velocity of propagation is assumed, then an 
estimate of the standoff distance can be made. Although this technique is 
impractical due to the extremely small time interval, it serves to illustrate the 
plasma density measurement concept. The velocity of energy propagation in the 
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plasma varies with the electron density profile according to (1.2.5). However, if 
the electron density variation is linear with distance, then the average group 
velocity is approximated well by the convenient value c/2, where c is the free 
space velocity of light. Consider a linear electron density profile and its 
corresponding dielectric profile as shown in figure 1 .2.2. 



e r = i- 




Z 


Figure 1 .2.2 Linear electron plasma density profile and the corresponding 
dielectric profile. 


The value of distance z p is the location of the plasma turning point. It is 

assumed here that only a single EM frequency (monochromatic) is launched 
into the plasma. The group delay t d within the plasma may be obtained 

analytically by integrating the inverse of the group velocity over distance. 



( 1 . 2 . 6 ) 


(1.2.7) 
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This shows that the virtual distance to the turning point is twice the actual 
distance. Therefore the average group velocity is c/2 in an electron plasma 
whose density profile is linear with distance. The CFD predicted profile shown in 
Fig. 1.2.1 is not exactly linear, but it is close enough to warrant the use of c/2 as a 
first approximation to the average velocity. 

It has been predicted by CFD that electron densities in the range of 10 12 - 
10 15 e'/cm 3 are expected to be encountered in the measurement range of 0-15 
cm from the heat tile during the aeropass. In order to measure this electron 
density range, frequencies of 9-284 GHz would be required. Due to the MRIS 
system space and weight constraints, as well as practical microwave and 
millimeter wave sources, four center frequencies have been selected by NASA: 
20, 44, 95, and 140 GHz. Thus, only four densities will be measured. 
Measurement of these four densities will be performed sequentially and 
repeated continuously for the first 500 seconds of the aeropass, and therefore 
can be tracked as a function of time. This information should provide enough 
data to accurately approximate the electron density profile as a function of time. 

In order to launch an EM wave into the plasma, an antenna is mounted 
behind the heat tiles of the AFE vehicle. The propagation path is illustrated in 
Fig. 1.2.3. The thermal protection heat tile has a dielectric constant e r = 1.36. 

Thus, a small dielectric discontinuity occurs at the vacuum-tile interface which 
will cause a small strength spurious reflection. However, there is a thin Reaction 
Cured Glass (RCG) coating on the outer surface of the heat tile whose dielectric 
constant e r is 4.8. Therefore, a significant spurious reflection will occur at the tile- 

RCG interface. This reflection is considered interference and represents a 
significant obstacle for the plasma density measurement. 
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Figure 1 .2.3 Microwave Reflectometer Ionization Sensor (MRIS) propagation 
path. 

The simulation results that will be presented for the two radar systems 
developed in this dissertation will be in the context of measuring electron plasma 
densities as an example of measuring target range in frequency dependent 
inhomogeneous media. However, the systems are much more general and 
application dependency of their utility is not to be inferred. 

1.3 Outline of Dissertation 


This dissertation is organized as follows. In Chapter 2 a method for 
modelling EM planewave propagation in inhomogeneous media is derived. The 
derivation is based on the theory of transmission lines. An inhomogeneous 
medium is modelled by a cascade of thin slabs where each individual slab is 
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homogeneous. The boundary value problem for the electric and magnetic fields 
is solved at each slab boundary. It is demonstrated that for an exponential taper, 
as the thickness of the slabs approaches zero, the error in the fields generated 
by the model asymptotically approaches zero. The model is then used to 
illustrate the behavior of EM radiation in electron plasmas. 

Chapter 3 presents the derivation of the frequency-stepped double 
sideband suppressed carrier (DSBSC) radar system. With certain modifications, 
this chapter appeared in [Ybarra2]. The newly derived system is an extension of 
the amplitude modulated (AM) CW technique [Nilssen] that was developed in the 
late 1 950's as a result of the stimulated interest in radar from World War II. 
However, the AM CW technique was abandoned because severely erroneous 
target range estimates were obtained when an interfering reflection was present. 
The main contribution of this chapter is the derivation of the DSBSC system 
which includes the effect of an interfering reflection. It is shown that the 
interference can actually be used as a reference allowing accurate target range 
estimates to be extracted from the radar data, independent of signal to 
interference ratio. Simulation results for the DSBSC radar system are then 
presented in the context of measuring electron plasma densities. 

The frequency-stepped continuous wave system is developed in Chapter 
4. Several versions of this technique are currently in use [Artech 87, Bird, 
Daviesl ,2,3, Wehner]. The basic system that is examined in this chapter is a 
generalized system formulation of the Hewlett-Packard HP-8510 network 
analyzer [HP]. The function of this chapter is to present a detailed description of 
the mechanism by which the CW radar data is gathered. The objective of the 
stepped CW system is to extract frequency response samples from the 
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propagation path. Once the data is obtained, signal processing must be applied 
to generate the target range estimate. The standard approach is to take the 
Inverse Fast Fourier Transform (IFFT) of the frequency domain data resulting in a 
time domain impulse response. Peaks in the impulse response correspond to 
reflections whose time delay corresponds to the target range. Unfortunately, the 
range resolution offered by the IFFT is inversely proportional to the total 
bandwidth of the frequency range spanned by the frequency steps. This 
resolution handicap of the IFFT approach is the motivation for Chapter 5 which 
presents the derivation of an optimal signal processing algorithm for the stepped 
CW data. The range resolution enhancement offered by this new algorithm is 
clearly demonstrated by applying it to two different sets of physical 
measurements. The first set of data was obtained by measuring the reflection 
coefficient at the input of a coaxial verification standard [HP] terminated in a 
matched load by an HP-8510 network analyzer. These measurements were 
performed in the ECE Dept, microwave laboratory at NCSU. The second set of 
measurements were obtained at NASA Langley Research Center (LaRC). An 
HP-8510 network analyzer was used as a microwave source driving a test 
apparatus for the MRIS instrument. This apparatus is described in detail in 
Chapter 5. The optimal signal processing algorithm is applied to the MRIS test 
data and the results illustrate significant resolution enhancement over IFFT 
processing. 

Finally, several future research objectives are proposed in Chapter 6 
along with a final statement of the conclusions that have been drawn from the 
research comprising the dissertation. 
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Chapter 2 


A Model for Electromagnetic Propagation in Inhomogeneous Media 

In this chapter a method for modeling plane wave propagation of 
electromagnetic (EM) radiation in inhomogeneous media is derived and 
illustrated. The model uses cascaded slabs of homogeneous media and is 
based on the theory of EM wave propagation in transmission lines. Several 
applications have utilized the concept of cascaded slabs to model 
inhomogeneous media. It is well Known that waveguides are dispersive 
[Collin], One way of minimizing phase distortion due to dispersion is to implant 
slabs of dielectric in the waveguide [Keuster]. In the case of a single dielectric 
waveguide, analytical solutions are possible only for a few specific permittivity 
profiles in simple geometries. In general, the problem can only be solved 
numerically. Kuester examines the numerical problem for arbitrary dielectric 
profiles. However, his numerical solution is subject to the constraints that all 
media have constant permeability p 0 , and that the relative dielectric constant 
satisfies £ r > 1 for all media. The model developed in this chapter is not bound 

by such constraints. One of the prevalent features of electron plasma is that its 
relative dielectric constant always satisfies e r ^ 1. Hence, a more general model 

is necessary for simulating propagation in plasma. In [Richmond] cascaded 
slabs of dielectric are used to model the effect of a temperature gradient on 
radomes. A very large temperature gradient in the radome of an aircraft often 
results from hypersonic flight through the atmosphere. The outer surface of the 
radome becomes hotter than the inside, resulting in a continuous variation in 
permittivity even if the radome was designed as a homogeneous structure. 
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Exact solutions in closed form are available only for a few special cases 
including linear and exponential inhomogeneities. In order to solve the radome 
problem, Richmond models the radome by a cascade of slabs, each of which is 
homogeneous. The procedure utilized is a step-by-step numerical integration 
of the first few terms of the power series expansion for the fields inside each 
slab. The technique is shown to produce results closely agreeing with 
experimental measurements. However, the case of frequency dependent 
inhomogeneities is not addressed and it is not clear how a reformulation could 
be made to incorporate frequency dependence. One of the primary 
contributions to this work from [Richmond] is the derivation of the relationship 
between slab thickness and modelling error. This topic will be discussed further 
in section 2.2. 

Since the model developed in this chapter is intended for simulating EM 
plane wave propagation in any inhomogeneous medium, it must be capable of 
incorporating frequency dependence in both permeability p and permittivity e. 
There are several methods currently in use for modelling plane waves in 
plasmas which incorporate frequency dependence in their non-uniform 
permittivity. The majority of these methods approach the problem in terms of 
finding the solution to a Fredholm integral equation of either the first or second 
kind. The unknown profile index (either e(x) or N e (x)) to be determined is 

contained in the integrand of one of the terms of the equation. Such 
approaches are termed profile inversion techniques [Ahn, Eden, Balanis1,2]. 
These profile inversion techniques appear to be as applicable to measuring 
inhomogeneous plasma densities as the cascaded slab model developed in 
this chapter. However, the solution of the Fredholm Integral equation for 
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arbitrary density profiles requires a linear approximation [Balanis2] as well as 
extensive computation. Therefore the need for developing other techniques for 
solving the plasma index profile problem exists. Reducing computational 
complexity while maintaining modelling accuracy is one of the primary 
motivations for current research in many engineering fields. The cascaded slab 
model developed in this chapter is implemented using the C programming 
language and incorporates the use of both recursive data structures and 
recursive programming techniques. The use of recursive programming 
techniques allows the simultaneous solution of all field quantities throughout 
the inhomogeneous medium. This leads to significantly reduced computation 
time over conventional programming techniques. 

Other techniques for plasma density determination such as those based 
on the WKB approximation [Ahn, Bahar, Daviesl ,2,3] are capable of providing 
accurate profile density estimates only for certain profile shapes. The 
cascaded slab model developed in this chapter requires no constraints on the 
shape of the profile index for accurate characterization. 
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2.1 Electromagnetic (EM) Characterization of Inhomogeneous 
Media 


Consider a vertically polarized EM plane wave propagating in a medium 
with permittivity e 1 as illustrated in Fig. 2.1.1. 



Figure 2.1 .1 . Plane Wave Propagation 

It is assumed here that the electric field intensity varies only in the spatial 
direction x and the magnetic field varies only in the y direction. This assumption 
may be expressed analytically by the relations 
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Assuming time harmonic fields (sinusoidal), the time reduced form of Maxwell's 
two curl equations (Harrington) is 


dE* 

dz 


= -jtOJUHy 


dH, 

dz 




(2.1.3) 
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These equations are similar to the differential equations governing the 
relationship between voltage and current for lossless transmission lines, 

f = S = - juCV (2.1.4) 

Hence, the theory used to describe propagation in transmission lines applies 
equally as well to plane-wave propagation. Now consider the discontinuity at 
the boundary between two different media each having different dielectric 
constants as shown in Fig. 2.1.2. This condition is similar to the discontinuity 
that arises when two different transmission lines with different characteristic 
impedances are connected to each other as shown in Fig. 2.1 .3. 



Figure 2.1 .2. Plane wave incident normally on a dielectric discontinuity. 
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Figure 2.1.3. Transmission line discontinuity 


The analogy between transmission lines and plane wave propagation (normal 
incidence) is summarized in Table 2.1.1. Note that the mismatch at the 
discontinuity is represented by a reflection coefficient for each case. 


Table 2.1.1. Plane Wave-Transmission Line Analogy [Brown]. 



Symbol or Equation 

Plane-Wave 

Quantity 

Symbol or Equation 

Voltage 

V 

Electric field Intensity 

Ex 

Current 

1 

Magnetic field intensity 

Jjy 

Inductance 
per unit length 

L 

Permeability 

p 


C 

Permittivity 

e 

Characteristic 

impedance 

yu 

n 

Intrinsic impedance 

EBiM 

Phase-shift constant 

P = covTC 

Phase-shift constant 

P = itifixe 

Velocity of 
Propagation 

v- 1 

Velocity of propagation 

V= 1 

VEC 


Reflection Coefficient 

r L = ZljlZo 

L Zl + Zo 

MfillSi 

1 


Incident wave power 

ihiihh 

11911 

Incident wave 
power density 
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Using the relationships provided by table 2.1.1, the problem of plane wave 
propagation through any number of discontinuities may be solved. Consider 
the system shown in figure 2.1.4 which is a cascade of multiple slabs each of 
which is homogeneous but whose combined structure is inhomogeneous. This 
problem may be solved by solving the equivalent transmission line problem 
shown also. Therefore, solving general transmission line networks is focussed 
upon first, and the problem of plane wave propagation is examined later in this 
section. 




Figure 2.1.4. Modeling propagation through multiple layers of media as 
transmission line sections. 


Consider the basic problem of simulating pulse transmission through a 
loaded transmission line. Assume that the pulse of interest is bandlimited with 
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a cutoff frequency of fc ■ The pulse response can be obtained by first computing 
the frequency response of the network and then multiplying the frequency 
spectrum of the pulse by the frequency response of the transmission line 
network. Finally, the inverse Fast Fourier Transform (IFFT) of this product yields 
the time domain pulse response. Therefore, as a first step in calculating the 
frequency response of the network, the network response to a single sinusoid of 
frequency fo is analyzed. Consider the loaded transmission line connected to 
the generator V$ through a source impedance Z$ as shown in figure 2.1 .5. 



x = 0 V(x), l(x) X = L 

Figure 2.1 .5. Generator connected to a loaded transmission line. 


The voltage and current at any point on the transmission line can be obtained 
from the following expressions (LePage] 
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i(x) = 


-yx 1 " r L e 


-2yO-x) 


z 0 + z s 


1 - r s r L e 


-2yL 


(2.1.6) 


In the above expressions the propagation constant y is given by 


Y= V(r + jcol)(g + jtoc) 


(2.1.7) 


and the characteristic impedance z Q of the transmission line is given by 


V g + jtoc 


(2.1.8) 


where r,l,c, and g are the distributed resistance, inductance, capacitance, and 
conductance of the transmission line respectively. The expressions for the load 
and source reflection coefficients are 


p _ z l~ z o 

L z l +z o (2.1.9) 

p _ z «~ z o 

s z s +z o (2.1.10) 


The expression for v(x) includes the superposition of all waves reflecting from 
the source and load mismatches. This can be seen by a Taylor series 
expansion of (2.1.5) 


v(x) = 1 e vx + r,e' Y(2L ' x) + r, r e V(2L+X) + 

Z- + Z L LS 


7 + Z 

A) ' 


2 2 


f l r s e Y(4L x) +r L r g e Y(4L+X) + ... ] 


( 2 . 1 . 11 ) 
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To obtain the shape of the pulse at the load, v(L) is evaluated at frequencies 
from f=0 to f=fc in discrete steps where f c is the cutoff frequency of the 
bandlimited pulse. Now consider the case where the boundary voltage and 
current are known on a section of transmission line as shown in figure 2.1 .6. 


v(0) 

1 ( 0 ) 



x=0 


v(x), i(x) 


x=L 


Figure 2.1 .6. Section of transmission line with boundary voltages and currents 


From (2.1 .5) it can be shown that knowledge of the voltage at x=0, v(0), allows 
the voltage at any point along the line v(x) to be computed from the relation 


v(x) 

v(0) 


-yx 1 + r L e" 2y(LrX) 

e -T77^ 


( 2 . 1 . 12 ) 


Similarly, the current at any point along the line i(x) may be computed by 


i(x) 

i(0) 


-YX 1 ■ ^L e 


■2y(l^x) 


1 - r L e 


-2 Y l 


(2.1.13) 
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Thus, using (2.1.12) and (2.1.13) the voltage and current can be evaluated at 
any point on the transmission line given the boundary voltage and current. With 
the above preliminaries, a methodology for solving the simple network shown in 
figure 2.1.7 is presented. In this figure, terminal planes, or nodes have been 
labeled ni through ns. To solve this network, that is to obtain the voltage and 
current at each node and at any location within the network, consider equation 
(2.1.5). This equation suggests that if the impedance at node ni were known, 
then the voltage and current at any location along the line can be calculated 
using knowledge of the source voltage and impedance. The first step is to 
obtain the impedance at ni. This impedance is seen to consist of the parallel 
combination of the impedances looking into ns and n 2 from n-|. 



Consider the impedance looking into a single loaded transmission line as 
shown in figure 2.1 .8. 
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Figure 2.1 .8. Input impedance of a loaded transmission line. 


This impedance can be obtained from the relation [LePage] 


Z in (x) = 


i + r^- 2 ^ 
1 - r i e' 2y(L ' x) 


(2.1.14) 


The first step is to calculate the impedances looking into nodes 03 and n4 from 
n2- The parallel combination forms the impedance at n2- The impedance at ni 
is calculated by the parallel combination of the impedances looking into 02 and 
ns. Therefore, the following methodology may be used to solve the network. In 
the first pass, starting from the three loaded end nodes, the impedances are 
calculated, and the parallel combination of these impedances at the parent 
node forms the parent node impedance. Working backward in this manner, the 
impedance at the root node (ni in the example) is calculated. Using (2.1.5) and 
(2.1.6) the voltage and current at the root node ni may be calculated. Then 
using (2.1.12) and (2.1.13) and the boundary voltages and currents calculated 
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at the parent node, the voltage and current at each node in the network can be 
calculated. 

In the special case of propagation through layers of different media, there 
are no branches. The transmission line analogy is illustrated in Figure 2.1.9. 



load 


Figure 2.1 .9. Network of cascaded sections. 

In this case the impedance at node ns is computed first. Next, the impedance at 
node n 4 is computed and so forth until the total impedance looking into the 
network is obtained. In the next phase, the voltages and currents are computed 
starting from the source, node ni, and then moving towards the load. 
Advanced computer programming techniques which utilize recursion can be 
implemented to efficiently perform these computations [Sobelman]. 
Transmission lines with a hundred thousand sections can be solved in only a 
few minutes of CPU time on a DEC 3100 workstation. It was shown by (2.1 .3) 
and (2.1.4) that the plane wave propagation problem is equivalent to an 
analogous transmission line problem. Utilizing the relationships between 
transmission line parameters and those for plane wave propagation as 
illustrated in table 2.1.1 allows the fields within a system of cascaded slabs to 
be computed. In the next section, the accuracy of modelling inhomogeneous 
media with cascaded slabs is analyzed. 
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2.2 Accuracy of Modelling Inhomogeneous Media with Cascaded 
Slabs 

In this section the accuracy of modelling an inhomogeneous medium 
with cascaded slabs is analyzed. Consider the transmission line illustrated in 
figure 2.2.1 which has an exponentially tapered characteristic impedance. 


Z Load 


L 

x=0 x=L 

Figure 2.2.1 Exponentially tapered transmission line. 

As a specific example, let the source have an impedance of 100 fl, the load an 
impedance of 500 Q, and total length of the exponentially tapered line be 10 m. 
The exponential taper can be created with two wire cable whose geometry is 
shown in figure 2.2.2. 
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D 

Figure 2.2.2 Parallel conductor transmission line geometry. 

The distance separating the pair of conductors is D and the radius of each 
conductor is d. The characteristic impedance of this transmission line is [Brown] 

Z 0 = 120 ln( 2D-) (2.2.1) 

For the continuous exponential taper the normalized impedance anywhere 
along the line Z(x) is related to the geometry and load impedance through the 
relation 

Z(x) =e 1 L l " Z| - oadl (2.2.2) 


The exact differential equation relating the reflection coefficient to location along 
the line is [Collin] 


dT = 
dx 


)2pr - 1 ( i- r 2 ) 


d In Z(x) 
dx 


(2.2.3) 


The exact solution to this Riccati differential equation for the input reflection 
coefficient is 


r i = 


A sin ( BL/2) 


B cos (BL/2) + j 2P sin ( BL/2) 


A = !^, B= V40 2 -A 2 


(2.2.4) 
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The normalized input impedance may be calculated from the input reflection 
coefficient using the relation 


i +r, 


z in = 


i -r 


in 


in 


(2.2.5) 


In order to test the model, sections of two wire transmission lines were 
cascaded. A comparison of the normalized input impedance as a function of 0 
obtained from the cascaded slab model and the exact solution to the Riccati 
equation (2.2.4 - 2.2.5) is shown in figures 2.2.3 and 2.2.4. The electrical length 
of the line 0 is related to the wavelength A and line length L through the relation 


0 = PL = 2il 
A 


( 2 . 2 . 6 ) 


The best comparison is shown in figure 2.2.4 which is an enlarged view of 
figure 2.2.3. The simulation performed with 500 sections could not be 
distinguished from the exact solution. The conclusion is that the cascaded slab 
model produces the exact solution asymptotically as the number of sections is 
increased. The reason for this asymptotic exactness is that the model includes 
the effects of all reflections, not a simple approximation where second and 
higher order reflections are neglected. 
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Figure 2.2.3 


Figure 2.2.4 


Exponential Taper 

Exact versus Computer Model (PropMod) 



Comparison of the normalized input impedance produced by the 
cascaded slab model and the exact solution to the Riccati 
differential equation. 


Exponential Taper 

Exact versus Computer Model (PropMod) 



Detailed comparison of the normalized input impedance 
produced by the cascaded slab model and the exact solution to 
the Riccati differential equation. 
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From this single evaluation of the accuracy of modelling an exponentially 
tapered inhomogeneity, it remains unclear as to how many slabs are necessary 
to model accurately an arbitrary inhomogeneous medium. It has been shown 
[Richmond] that the error in modelling a frequency independent 
inhomogeneous medium with cascaded slabs is bounded by the two 
expressions 

e(h) < khV^ (2.2.7) 

« ■ ^i 1/4 

e(h) <kh|e r Ve r | (2.2.8) 

i 

where h is the slab thickness, k is the measurement uncertainty, fy is the 

derivative of permittivity with respect to distance, and e(h) is the modelling error 
relative to the measurement error. The error bound (2.2.7) simply shows that 
reducing the step size reduces the modelling error. The more important of the 
two equations is (2.2.8) which shows that the step size h should be reduced 
when the slope of the permittivity profile is large (steep). This result obeys 
intuition since any spurious reflection will lead to modelling error. Spurious 
reflections will occur when the model produces discontinuities that are not 
present in the medium being modelled. When the profile slope is large, more 
slabs are needed to provide a smooth transition from one slab to the next. The 
derivation provided by Richmond considers only media with permittivity greater 
than one. Hence, the equations bounding the modelling error (2.2.7) and 
(2.2.8) cannot be applied to plasma profiles, where the permittivity is less than 
one. No error bound appears to exist in the literature for using cascaded slabs 
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to model inhomogeneous plasmas. The derivation of such a bound remains as 
future research stimulated by the work leading to this dissertation. 

In order to model the plasma, the electron density profile is sliced into 
slabs of varying thickness, each characterized by its complex permittivity 
determined by (1.2.3). A dynamic slab thickness allocation algorithm should 
allocate slab thickness inversely proportional to the magnitude of the derivative 
of the dielectric profile with respect to distance as suggested by Richmond. This 
constraint is required to maintain smoothness in the dielectric profile of the 
model since the true inhomogeneous plasma has a continuous dielectric 
profile, and any sudden discontinuity in the dielectric will cause spurious 
reflections. With a plasma model, another consideration must be made. It is 
important that the modeling error be relatively small in the vicinity of the turning 
point. Therefore, considering slope and the turning point criterion, a dynamic 
slab width allocation algorithm was developed which considers both criteria in 
the allocation process. The algorithm is quite simple and is illustrated in figure 
2.2.5. The number of sections per unit distance decreases linearly with 
increasing permittivity for linear profiles. Only 85% of the total allocation of 
sections are shown in the diagram. The remaining sections are evenly 
distributed between the turning point and the end of the profile which is located 
at a standoff distance of approximately 1 5 cm (see figure 1 .2.1). 
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% of total sections 



Figure 2.2.5 Dynamic allocation of plasma model slab widths. 


Figure 2.2.6 shows what happens with this distribution for three differently 
shaped dielectric profiles. In the linear case, the slope is constant and 
therefore there is a linear increase in sections per distance as the turning point 
is approached. In the second case, the slope is near zero for nearly half the 
distance to the turning point and then the slope becomes increasingly larger 
(negative) causing the sections to compress nearer to the turning point. This is 
precisely what is desired, since the slope is steepest near the critical distance 
d c . The third case shows the resulting trade-off between slope and dielectric 

value. The slope is large near the front edge of the profile and approaches zero 
near the critical point. Thus, the algorithm allocates relatively few sections in 
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the first few permittivity bins. However, this results in thinner slabs than would 
occur for a linear taper due to the large, negative slope in the permittivity profile. 
Then as the critical point is approached, the number of sections is increased 
because the permittivity is decreasing. However, because the slope 
approaches zero with increasing distance, the slab width actually increases 
slightly as the plasma turning point is approached, and a slope/turning point 
trade-off occurs as desired. 

The total number of sections used in the slab model for accurate 
representation of the plasma was determined as follows. For the particular CFD 
profile used in the simulations, it was found that using more than 1000 sections 
made very little difference in any of the measurements. In order to assure an 
accurate plasma model for all measurement cases, 1500 sections were used in 
all simulations. Using 1500 sections resulted in section sizes ranging from a 
few microns to a few millimeters. 

Although the cascaded slab model is capable of modelling accurately 
planewave propagation in plasmas, it must be pointed out that there is one 
significant drawback. In a physical measurement environment where the 
transmitted EM energy is emitted near the plasma, such as in the MRIS 
experiment, the resulting EM energy is not composed of planewaves. The 
energy spreads with range and no account of this energy spreading is currently 
provided by the cascaded slab model. Hence, one of the future research 
objectives is to include the effects of spherical wave transmission through the 
layers of media. 
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Figure 2.2.6 Illustration of the algorithm which dynamically allocates slab 
thicknesses in the plasma model for various permittivity profiles. 
Tick spacing indicates the relative number of sections used in the 
plasma model. 
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2.3 Obtaining Impulse and Frequency Responses of Channels 
Composed of Reflective Targets in Inhomogeneous Media 


This section investigates the use of the cascaded slab model to obtain 
impulse and frequency responses of channels in which a reflective target is 
embedded in a dispersive, inhomogeneous medium. In each example, an 
electron plasma is used as the inhomogeneous medium. The electron plasma 
quantities used in the cascaded slab model are presented in table 2.3.1 [Chen]. 


Table 2.3.1 . Electron plasma relationships used in the cascaded slab model. 


Permittivity, n = electrons/cm 3 

e-eoEr-eof 1 n ) 
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As a first example, consider figure 2.3.1 which shows a constant density 
electron plasma sandwiched between two layers of air. A metal plate is placed 
against the second air layer. In this example, we are interested in the impulse 
response at the boundary of the first slab where it is assumed that a packet of 
plane waves is incident on the system. The electron density of the plasma is 
44.5*1 0 1 2 e'/cm 3 corresponding to a critical frequency of 60 GHz. 


AIR AIR metal 

e r =i PLASMA e r =1 P^te 



Figure 2.3.1 Channel composed of a constant density N e = 44.5*1 0 12 e'/cm 3 
electron plasma sandwiched between two slabs of air and 
terminated in a metal plate. 

Two cases are now presented. In the first case, the impulse response 
over a bandwidth of 20 GHz at a starting frequency of 40 GHz is considered. For 
the second case, the same bandwidth is used but the beginning frequency is 
increased to 100 GHz. In both cases 512 frequency samples were obtained, 
which provides 1024 points in the resulting impulse response computed from a 


36 



512 point IFFT. The impulse response corresponding to the first case is shown 
in figure 2.3.2. The time domain data has been converted to the distance 
domain using the free-space velocity of propagation. In addition, consideration 
of the two-way travel time has been provided in the distance scale. Since the 
frequency span is 40 to 60 GHz, a strong reflection is observed from the plasma 
whose critical frequency is 60 GHz. The range to the reflection is 20 cm, 
precisely as expected from the channel geometry (figure 2.3.1). 


START FREQUENCY 40 GHz, BANDWIDTH 20 GHz 


impulse response 



Figure 2.3.2 Impulse response of electron plasma sandwiched between air 
slabs and terminated in a metal plate. The starting frequency is 
40 GHz. 
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In the second case, the frequency span is from 100 GHz to 120 GHz. The 
resulting impulse response is shown in figure 2.3.3. 


START FREQUENCY 100 GHz, BANDWIDTH 20 GHz 


impulse response 



Figure 2.3.3 Impulse response of electron plasma sandwiched between air 
slabs and terminated in a metal plate. The starting frequency is 
100 GHz. 

Most of the energy passes through the plasma because its critical 
frequency is 60 GHz. However, there are minor reflections due to the dielectric 

discontinuities at the air/plasma interfaces. At 100 GHz, the dielectric constant 
of the plasma is e r = 0.64. The resulting reflections are shown in figure 2.3.3. 

The distance to the main reflection is shown to be 40 cm which corresponds to 
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the distance to the metal plate as indicated in figure 2.3.1. The velocity of 
propagation within the plasma has been taken into account. 

In the next example, an analysis of propagation through the CFD electron 
density profile is presented. Consider the CFD profile shown in figure 2.3.4, 
which has been scaled so that the critical electron density at 20 GHz is located 
at a standoff distance of 4 cm. Note that the vertical axis is critical frequency 
and not electron density. The relationship between these two quantities was 
given in (1.2.4). 

Profile for 4 cm Standoff Distance at 
20 GHz Critical Frequency 



Figure 2.3.4 Scaled CFD electron density profile for 4 cm standoff at 20 GHz 
critical frequency. 

This profile is placed outside the heat tile system of the AFE spacecraft which 
was shown in figure 1 .2.3 so that the standoff distance is zero at the outer 
surface of the tile. The frequency response of the heat tile system and CFD 
profile is shown in figure 2.3.5. The reflected energy is added to the transmitted 
wave and forms an interference pattern. The sudden drop in reflected energy at 
30 GHz corresponds to the maximum critical frequency in the CFD profile. The 


39 



periodic behavior of the frequency response beyond 30 GHz is due to the 
presence of the RCG only and its period is inversely proportional to the distance 
from the antenna aperture plane to this source of reflection. 


Plasma plus Tile/RCG Frequency Response 
at the Source 



Figure 2.3.5 Tile/RCG and plasma frequency response as seen from the 
MRIS antenna flange. 

Figure 2.3.6 shows a close-up view of the frequency response between 20 and 
30 GHz. The distance between the nulls in the frequency response decreases 
as the frequency increases in an "accordion-like" effect. This is due to the fact 
that the turning point is getting farther away from the source. 
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Plasma plus Tile/RCG Frequency Response Magnitude 

. Measured at the Source 

frequency response 24 GHz < f < 28 GHz 

magnitude 



Figure 2.3.6 Close-up of plasma frequency response showing reflections 
from changing critical densities. 

The next set of results provides great insight into the behavior of the 
plasma. Consider figure 2.3.7 which shows the impulse responses obtained 
using the cascaded slab model for the Tile/RCG CFD profile system using 64 
frequency steps at 64 MHz intervals. This corresponds to the proposed 
measurement sequence for the MRIS instrument. Each impulse response 
corresponds to a different starting frequency as indicated. The frequencies are 
stepped in 64 MHz increments up to 64 steps. Therefore, the bandwidth of each 
measurement sequence is 4.096 GHz (64 x 64 MHz). 
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Impulse Responses as the Starting Frequency 
is Changed (BW = 64x64 MHz) 



Figure 2.3.7 Impulse responses of the CFD profile as the starting frequency is 
changed. 

The reflection from the plasma turning point moves outward, following the 
profile, as the starting frequency is increased. The turning point reflection may 
exhibit 180° phase shifts as in the case for 24 GHz. At 28 GHz, the top of the 
profile is skimmed (see figure 2.3.4). At 30 GHz and beyond there are no more 
plasma reflections. The only reflection is from the RCG, and all impulse 
responses are basically the same. The impulse responses up to 40 GHz in 2 
GHz steps were computed and are shown in the graph. 

An interesting view of the data from figure 2.3.7 is provided by the three 
dimensional plot shown in Figure 2.3.8. In this figure, all 15 impulse responses 
are plotted up to 18 cm standoff as indicated. The reflections below 30 GHz 
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start out close to the tile outer surface and then move out and eventually 
disappear. As the frequency increases beyond 30 GHz, the impulse responses 
all flatten out with no reflections other than those resulting from the tile. 



18 40 


Figure 2.3.8 3 dimensional plot of impulse responses of Tile/RCG and CFD 

profile system as the starting frequency is changed. 


2.4 Summary 

A model for electromagnetic planewave propagation in inhomogeneous 
media was derived in this chapter and its performance illustrated in the context 
of modelling EM propagation in non-uniform electron plasmas. This model is 
based on the theory of transmission lines and represents an inhomogeneous 
medium by thin cascaded slabs, each of which is homogeneous but whose 
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combined structure is inhomogeneous. It has been shown that as the thickness 
of each slab is reduced, the response from the model asymptotically 
approaches the true physical response. The model can be implemented very 
efficiently using advanced computer programming techniques such as recursive 
programming which allow the use of hundreds of thousands of slabs while 
requiring little computation time. 

Several simulation results have been presented in the context of 
measuring non-uniform electron plasma densities that develop near the heat 
tiles of a space re-entry vehicle. Unfortunately, an appropriate plasma profile 
cannot be produced for ground testing. Hence, with the exception of performing 
an actual atmospheric aeropass, their measurement can only be performed by 
computer simulation. For this reason, the inhomogeneous media model 
derived in this chapter, which can simulate accurately the propagation of 
planewaves in electron plasmas, is especially significant. 
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Chapter 3 


The Frequency-Stepped Double Sideband Suppressed Carrier 

Radar System t 

3.1 Introduction to the Frequency-Stepped DSBSC System 

The double sideband (DSB) radar system derived and demonstrated in 
this chapter is a reformulation of the Amplitude Modulated CW radar system 
developed in 1962 [Nilssen]. Both systems extract target range information by 
observing the phase delay of the echo signal. The idea of using phase 
information to extract target range was first introduced in 1947 [Riden] during a 
time in which much effort was spent by scientists from several fields who joined 
in an effort to promote the development of practical radar systems. The AM CW 
radar technique has received little attention since its introduction in 1962, 
because the system requires the echo from the primary target to be much 
stronger than any other echo. The presence of an interfering reflection 
rendered the target range estimate useless. The DSB technique derived in 
section 3.2 can operate in the presence of an interfering reflection. It is shown 
that accurate target range may be extracted for virtually any signal-to- 
interference ratio (SIR). This immunity from weak target echoes is achieved by 
frequency-stepping the suppressed carrier and obtaining a sequence of phase 
measurements. The original formulation used only a single carrier frequency 
and for this reason was unable to mitigate the unwanted affects of interference. 

A field closely related to radar is remote sensing of the atmosphere. 
There exists a region of the atmosphere called the ionosphere which is 

t This chapter appeared, with certain modifications, in [Ybarra2j. 
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characterized by a large number of oxygen and nitrogen ions resulting from 
solar radiation. The density of these ions is great enough to affect the 
propagation of radiowaves. The ionosphere begins at an approximate altitude 
of 50 km and reaches a maximum ion density at approximately 80 km. It has 
several unusual EM wave propagation properties. One of these properties is 
that it reflects radiowaves with certain frequencies somewhat similarly to a 
conductor. This reflection property allows certain radio frequencies to "skip" off 
the ionosphere and travel literally around the globe. Because the ionosphere 
plays such an important role in radiowave communication, remote sensing of 
the atmosphere has remained a highly active area of research for more than 30 
years [Daviesl.2,3]. The primary process currently in use for measuring the ion 
density of the ionosphere is the ionosonde. lonosondes are performed from 
both terrestrial sites as well as from satellites. There are two primary methods 
of measurement: (i) measurement of the time delay of a pulse and (ii) 

determination of the stationary phase. The direct measurement of pulse delay 
is done by starting a timing device (e.g. cathode ray oscilloscope time base or a 
digital clock) when the transmitter is triggered and noting the time of the leading 
edge of the echo pulse(s) [Davies3]. Uncertainties arise in identifying which 
feature on the pulse envelope to use as reference. Even with a relatively 
square emitted pulse, ionospheric dispersion tends to produce a triangular or 
trapezoidal shaped echo. In some cases increased accuracy is obtained by 
using Gaussian-shaped pulses [Devlin]. However, time accuracy is limited to 1 
ps which corresponds to 0.3 km of virtual height. Because of this inherent 
accuracy limitation, the time delay method is impractical for short range radar 
systems. 
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The phase approach is based on the relationship between phase and 
time delay which is also the basis for the derivation of the stepped DSB 
technique presented in the next section. In this ionosonde method, a single 
frequency sinusoidal signal is transmitted into the atmosphere and the phase of 
the echo is compared to the phase of the sinusoid produced by the generator. 
From knowledge of the frequency and velocity of propagation, the distance to 
the reflection point can be estimated. Although this technique is impractical for 
target ranging in general radar applications for the very same reasons the AM 
CW radar failed, it provides another example of using phase measurements to 
obtain the range to a reflective target. 

There are two currently existing short range radar systems which are 
somewhat similar to the DSBSC technique and are used for automotive radar 
[Boyer, Grimes]: the Diplex Doppler method and the Sinusoidal FM method. In 
the Diplex Doppler method, the transmitter is time-shared between two closely 
spaced frequencies. The phase difference between the two slightly displaced 
Dopplers at the receiver depends on the range to the target. However, since 
the range is extracted from the phase difference in the Doppler signal, no 

information is obtained without significant target velocity. 

The Sinusoidal FM technique extracts the target range from the relative 

phase of the received signal which does not depend on target speed. 
Unfortunately, the Sinusoidal FM system is relatively complex to implement 
[Grimes]. The advantage of the DSBSC technique is that it is simple to 
implement, requiring minimal hardware and space, while providing accurate 
target distance estimates. An important additional advantage of the frequency- 
stepped DSBSC system over other techniques such as a stepped-carrier 
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technique that generates an impulse response and uses range-gating to 
separate the target from an interferer is its immunity to Doppler shift due to 
target motion. The relative Doppler immunity of the frequency-stepped DSBSC 
technique is due to the fact that the effective Doppler shift is at twice the 
modulation frequency which is much smaller than the microwave carrier. 

3.2 Derivation of the Frequency-Stepped DSBSC Technique 

The DSBSC distance measurement technique is based on the 
relationship between the phase response and time delay of a linear system. 
For a linear system, the time delay is the negative slope of the system phase 
response and is therefore, in general, a function of frequency. 

t - d6(cc) 

dto (3.2.1) 

In the time domain, the simplest realistic model for the impulse response of a 
propagation channel known to contain an interfering reflector and a highly 
reflective target is given by 

h(t) = Ap5(t-tdp) + Aj6(t-tdi) (3.2.2) 

where 5(t) is the Dirac Oeita function, and. *. dp and. f dj a*a the round trip time 

delays to the primary target and interfering reflector respectively. The constants 
A p and A f represent the combination of path loss and reflection coefficient 

associated with each reflector. This model assumes that a transmitted 
microwave energy burst will return primarily as two temporally concentrated 
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pulses. Simulation results presented later show that the result of using this 
assumption is accurate for the plasma measurement application. 

The frequency response associated with a channel whose impulse 
response is given by (3.2.2) is 

H(f) = A p e*i wt dP + Aje-j 10 ^ (3 2 3) 


where co = 27tf. Separating H(f) into its real and imaginary parts yields 

H(f) = Apcos(cctdp) + Ajcos(cotpi) - j [A p sin(wt<jp) + AiSin(cotdi)] (3.2.4) 


The phase of this frequency response function is 

9(f) = Tan j - - ( Apsin(cotdp) + 

\ A p cos(cot(jp) + AiCOS(wt(ji) J 


(3.2.5a) 


The phase response may be expressed in the alternate form 


0(f) = -27rftd P + tan' 1 


A|Sin(27rf(tdp - tpj)) 

Ap + A|C0s(27rf(t C jp ■ tpj) j 


(3.2.5b) 


The negative derivative of the phase response (3.2.5a) with respect to angular 
frequency w is the channel time delay and is given by 


^ / ApSin(cotdp) + AjSin(totdi) \ 2 

(Apcosfcotdp) + Ajcos(cot(jj)/ (3.2.6) 
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where 


P- 


C + D 


[ApCOS(<ot dp ) + AjCOS(cot dj )] (3.2.7) 

C = [ApCOS(wt d p)+AjCOs((ot d j)j[Apt d pCOS(cot d p)+ Ajt d jCOS(wt dj )] ^.2.8) 
D = [ApSin(tot d p)+AjSin(cct d j)][Apt d pSin(wt d p)+Ajt d| sin(wt d j)] ( 329 ) 


Using the definition 


a = 


Ap 


( 3 . 2 . 10 ) 


which is the signal-to-interference ratio (SIR), and elementary trigonometric 
identities, the expression for the channel time delay (3.2.6) may be simplified to 

_ CX 2 tdp+ tdi+Ct(tflp+ tdi)COS(27Tf(tdp- t<ji)) 

1+a 2 +2aCOS(27rf(V t <Ji)) (3.2.11) 


This function is periodic in frequency with period 


tdp-tdi (3.2.12) 

Thus, one mechanism for determining the distance to the primary target is to 

estimate the period of the group delay measurement. Since the time delay to 
the interfering RCG coating t dj is known a priori, the period (3.2.12) of the 

measured group delay yields the delay time to the plasma turning point t dp . 

Using an estimate of the average energy propagation velocity for the plasma 
enables t dp to be converted into a distance measurement. It will be shown 
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through simulation that this method of extracting the turning point location may 
lead to accurate range estimates. 

A previously reported method of extracting the turning point standoff 
distance using the frequency-stepped DSBSC approach [Hearn] is to average 
the phase measurements, or equivalently the time delay measurements. The 
average time delay may be derived analytically by integrating the time delay 
expression (3.2.6) over one period and dividing the result by the period. The 
expression for the time delay (3.2.6) has the form 

f( x ) - a+bcos(x) 

c+dcos(x) (3.2.13) 


where the following substitutions have been made: 

a = a 2 td P + tdi, b = a(to p + y, c = 1+a 2 , d = 2a, x=27rf(t dp - tdi)- 


Given the fact that c 2 £ d 2 for all a, the integral of f(x) over one period is 


r a+bcos(x) dx = 
I c+dcos(x) 


2(a - ^) (Vc2-d 2 tang) 

bx + d_ Tan -1 m 

d Vc2T^2 ' 


c + d 


7T 


-71 


(3.2.14) 


The proof of this antiderivative is presented in the appendix. Evaluation of 
(3.2.14) and dividing by the period 2 tt yields the average value f gv given by 


iav 




a- be. 

d_ 

Vc 2 - d 2 


(3.2.15) 
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Now resubstituting the coefficients a, b, c, and d yields the average time delay 
'dav 9 iven b V 


■>. (tdp+tdl)! 1 +a2 ) 

. _tdp + ldl , ” 1,111 ta '~ 2 

^ 2 ~ |a 2 -l| 


(3.2.16) 


This expression may now be simplified to yield the result 

td av = gjtdp + ^ + j^jT ^ 


(3.2.17) 


which may be further simplified to obtain the following final result 


^av ~ *dp ° > ^ ’ *dav *di a < ^ 


(3.2.18a) 


The result for the special case of a=1 may be obtained by substituting this value 
directly into (3.2.1 1). After simplification, the average time delay for this case is 


. _ tdp + tdi 

KteV n 


(3.2.18b) 


The analytical results derived here for the time delay obtained using the 
DSBSC technique may be interpreted as follows. If the plasma reflection is 
stronger than the reflection from the interfering RCG (a > 1), and if an accurate 
technique for extracting the average phase is employed, then the turning point 
may be established. If the SIR is less than one (a < 1), and phase averaging is 
employed, an erroneous estimate of the turning point will result. It has been 
analytically shown, however, that an accurate estimation of the turning point 
may be extracted from the phase measurements using the relation between the 
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group delay period (3.2.12) and the time delay difference t dp - t dj , regardless of 

the signal-to-interference ratio. This result shows the potential immunity of the 
DSBSC technique to a weak primary target reflection. 

The theoretical phase measurement, assuming two scatterers, has been 
shown to be 

0(f) = Tan^f " (MHfrfrp* + A iSin(cotdi)) \ 

\ Apcos(cotdp) + Ajcos(wtdi) ) (3.2.5a) 

Plots of the theoretical phase (3.2.5a) over a 4 GHz bandwidth, beginning at 18 
GHz, for various SIR'S are shown in figures 3.2.1 a-d for the two scatterer 
idealization of the propagation path of figure 1.2.2 with the plasma critical 
density at 9 cm for 20 GHz (an average group velocity of c/2 has been assumed 
for the plasma). 
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Figure 3.2.1 


phase (radians) 



(a) 


phase (radians) 



(b) 

Theoretical frequency-stepped DSBSC phase measurement 
assuming two scatterers and SIR >1. (a) 4 GHz bandwidth 
(b) 600 MHz bandwidth. 
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Figure 3.2.1 


phase (radians) 



(c) 


phase (radians) 



(d) 

Theoretical frequency-stepped DSBSC phase measurement 
assuming two scatterers and SIR < 1 . (c) 4 GHz bandwidth (d) 
800 MHz bandwidth. 
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Several important observations can be made from these graphs. First of all, 
figures 3.2.1 a-b show that if the interferer is absent, the SIR a=°°, and the phase 
is linear due purely to the primary target reflection. As the SIR deteriorates 
approaching unity from above, an oscillation about this linear phase occurs with 
increasing amplitude and may be seen most clearly in figure 3.2.1b which is an 
enlarged view of figure 3.2.1a. However, the periodicity of this oscillation does 
not change with SIR except in the extreme case a=°°. Figures 3.2.1 c-d illustrate 
a similar phenomenon for cases of SIR’S between 0 and 1. When a=0, the 
phase becomes linear once again, and this is due purely to the interfering RCG. 
Note that the periodicity of the oscillations, once again, does not vary with SIR. 
When the SIR is identically equal to 1 , the phase function will converge to a 
linear phase with an intermediate slope exactly half way between the extreme 
cases of a=0 and a=1 , exactly as predicted by the theory developed earlier. It 
has been shown analytically that the periodicity of the oscillation about the 
linear case has a period which is the inverse of the difference in time delays to 
each scatterer. This relation is indicated in figures 3.2.1a and 3.2.1c, and 
special note should be made that this period is constant, independent of SIR. 

The time delay may be calculated from the phase response by 
differentiating with respect to to and then may be converted into a distance 
measurement as a function of frequency using an average plasma velocity. 
Using an average plasma velocity of c/2, several theoretical distance 
measurements using the same 4 GHz bandwidth as was used in figure 3.2.1 
are shown in figure 3.2.2 for various SIR'S. 
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distance (cm) 
20 



frequency 
22 (GHz) 


Figure 3.2.2 Theoretical distance measurements as a function of frequency 
for the frequency-stepped DSBSC system. 


It has been proven that if the SIR is greater than one, then the average value of 
the distance function for any SIR is the primary target distance. Furthermore, it 
has been proven that if the SIR is less than one, then the distance function has 
an average value equal to the distance to the interfering RCG. Averages were 
taken for each of the distance functions shown in figure (3.2.2). When the SIR is 
such that the deviation from the mean value is small (a near °° or a near 0), 
then only a few points need to be averaged in order to obtain an accurate 
estimate of the average. However, when the SIR approaches unity from above 
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or below, many points are required for an accurate estimate of the average. 
This is due to the near cusps that develop at the function extrema. The 
periodicity is shown to contain the standoff distance information as expected 
from the theoretical derivation. An important point must be made regarding the 
process of averaging the distance measurements. Accurate estimates require 
that an integer number of half periods are averaged or that a large number of 
half periods be averaged. When the target and interferer are far apart, several 
periods are present. When the target is near the interferer, the period 
approaches zero which may cause an inaccurate distance measurement using 
the distance averaging approach especially if the SIR is near unity and only a 
relatively small number of measurements are averaged. In this case, it is 
fortunate that the distance information may still be extracted from the periodicity 
of the phase measurements, even when the oscillation is weak. 

3.3 The Frequency-Stepped DSBSC Radar System Block Diagram 
and Principles of Operation 

The Frequency-Stepped DSBSC system simultaneously transmits two 
sinusoids centered around the suppressed carrier. The differential phase shift 
between these two sidebands induced by the propagation channel is then 
measured at the receiver producing an incremental estimate of the slope of the 
phase response. This estimate of the phase response slope is the estimated 
group delay. A measurement is made at N discrete frequencies, spanning a 
total bandwidth of NAf, where Af is the frequency step size. An implementation 
of the DSBSC system is shown in figure 3.3.1. 
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Figure 3.3.1 Block diagram of the frequency-stepped DSBSC radar system. 


59 




The transmitter consists of an upper single sideband frequency-stepped carrier 
generator whose output is modulated by another sinusoid of frequency f M . The 

received echo is passed through a square-law device and bandpass filtered to 
extract the component at frequency 2f M produced by the square-law device 

since it is this tone that contains the phase differential 0 2 -0 1 where 0 2 and 0 1 

are the channel induced phase shifts of the upper and lower sidebands 

respectively. This signal is then simultaneously mixed with an in-phase and 
quadrature sinusoid at frequency 2f M , phase coherent with the modulator in the 

transmitter. After low-pass filtering, the result is a pair of quadrature signals 
from which the phase differential 0 2 -0 1 may be extracted via a four quadrant 
arctangent operation. The round-trip group delay estimate t d is obtained from 

the phase difference by 


» _ _ d0(<o) a 02 - 9 i _ 82-61 
J die W 2 -U 1 4jrfM (3.3.1) 

The maximum unambiguous range corresponds to a 2n variation in 0 2 -0 1 and 
is given by 


FW= 2 Vflta= it 


(3.3.2) 


using an average plasma velocity of c/2. Since it is desired to measure electron 
densities with standoff distances from 0-15 cm, we chose to use f M = 125 MHz in 

our simulations which provides a maximum unambiguous range of 
approximately 30 cm from the antenna aperture. The range equation (3.3.2) 
dictates an upper bound on f M . The smaller the value of f M , the better the 
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differential phase approximates the true slope. However, as the value of f M is 

reduced, the phase differential corresponding to a particular standoff distance 

decreases. Therefore it becomes increasingly difficult to resolve small standoff 
distances as f M is decreased. 

The block diagram of the DSBSC system presented in figure 3.3.1 shows 
the system exactly as implemented using CAPSIM, a hierarchical block diagram 
communication and signal processing simulation environment. The block 
stairgen creates a sequence of samples with a staircase amplitude which drives 
the digitally-controlled oscillator DCO to produce the discrete frequency sweep 
of N steps over a bandwidth of NAf. The block in figure 3.3.1 labelled Target in 
inhomogeneous medium utilizes the inhomogeneous media model derived in 
Chapter 2. 

It still remains to be shown whether the two scatterer theory developed in 
section 3.1 is applicable to the much more complex vacuum/tile/RCG/plasma 
propagation path which was shown in figure 1.2.2 (pg. 8). Simulation results 
will now be presented which show that the theory is indeed applicable. 

3.4 Simulation Results for the Frequency-Stepped DSBSC System 


In order to illustrate the performance of the frequency-stepped DSBSC 
system, consider the CFD electron density profile shown in figure 1.2.1 (pg. 5) 
scaled so that the critical density at 140 GHz is located at a standoff distance of 
9 cm. Using N=64 measurements with a frequency step-size of Af=64 MHz, the 
DSBSC system produced the distance measurements shown in figure 3.4.1 . 
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Figure 3.4.1 Distance measurements using the frequency-stepped DSBSC 
radar system with and without loss caused by plasma collisions. 

When no loss due to collision is included, the group delay is captured by the 
primary target because the reflection from the plasma is stronger than the 
reflection from the RCG interferer. However, when loss is included in the 
simulation, the reflection from the RCG is much stronger than the plasma 
reflection and the group delay measurement is captured by the interfering RCG. 
The SIR in this case is approximately 0.1 (-20 dB). This SIR value was obtained 
by using a stepped CW technique and comparing the relative strengths of the 
plasma and interfering RCG reflections in the time domain. The averaging 
technique described in [Hearn] would work if no collision were present, 
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although producing an error of 0.5 cm (i.e. the average value of the no-collision 
curve is 12.70 cm). However, the averaging technique will fail in the more 
realistic case which includes path loss due to collision. Fortunately, the 
information containing the location of the plasma target is present in the data 
periodicity whether or not loss is present. It was shown in (3.2.12) that the 
fundamental frequency of the measurements shown in figure 3.4.1 is the 
difference in the time delays to the primary target and interferer assuming the 
simple two scatterer model. There exist several ways to extract the fundamental 
frequency from the N data points available from a given measurement sweep. 
The next section examines the use of autoregressive analysis on the group 
delay measurements to obtain the target range. 


3.5 Autoregressive Signal Processing of the Frequency-Stepped 
DSBSC Data 

One approach to extracting the fundamental frequency from the N data 
points available from a given measurement sweep would be to take an N-point 
Discrete Fourier Transform (DFT) of the data. The time resolution resulting from 
an N point DFT of the frequency-stepped DSBSC distance data is 

At = --- — 

2NAf (3.5.1) 


Using an approximate average plasma group velocity of c/2, the one-way 
distance resolution is 


= H At = 

2 2 8NAf (3.5.2) 
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Using N = 64 and Af = 64 MHz, a distance resolution of Ad = 0.9155 cm results. 

An alternative analysis technique is now presented for estimating the 
standoff distance from the periodicity of the distance measurements. It is based 
on autoregressive time series analysis [Kay81 ,88] which has been applied to 
high resolution spectral estimation. As will be shown, time series analysis may 
provide a significant resolution improvement in standoff distance estimates over 
that obtainable using the Fourier approach. 

Consider the near periodic sequence of distance measurements shown 
in figure 3.4.1 to be approximated by the difference equation 

3 

d(n) = £ akd(n-k) + u(n) n = mAf 0<m<63 

k=i (3.5.3) 

where d(n) is the sequence of distance measurements, a k are the 
autoregressive coefficients, and u(n) is the innovation or input driving function. 
Thus, d(n) is being modeled as the response of a linear system whose output is 
the innovation plus a linear combination of the past 3 outputs. Such a model is 
called an autoregressive (AR) model of order 3. A third order autoregressive 
model was chosen because a single peak in the time response is expected 
which requires two complex poles in the model. The other parameter is used to 
model deviation in the phase measurements from a pure sinusoid. In 
conventional time series analysis d(n) would be, typically, discrete samples of a 
continuous-time function d(t). Here, the d(n) are discrete samples of a 
continuous-frequency function d(f). If the discrete-time Fourier Transform of 
(3.5.3) is taken, a continuous complex function of time results. 
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(3.5.4) 


3 

□(gjAiotam) = £ akD(ei Awtdm )e‘l kAtll3tdW + U(ej Atotdm ) Aw = 2nAf 

k=1 


where t djff = t dp - t dj . The system time response (output D over input U) may be 
obtained from (3.5.4) resulting in 


R(W = 


1 - 


1 


3 


I 


a k e-ikAwu* 


(3.5.5) 


which is once again a complex function of time. In order to illustrate the value of 

this approach, the time series model (3.5.3) was applied to the collision 
distance data of figure 3.4.1 . The autoregressive coefficients a k were extracted 

using a conventional least-squares approach [Kay81,88]. These coefficients 
were then substituted into (3.5.5) and the resulting squared magnitude is shown 
in figure 3.5.1. The delay estimate is 1.203 ns as shown in the figure and 
translates into a standoff estimate of 9.0225 cm, producing an error of only 225 
pm. It should be noted that (3.5.5) is a continuous function of time and does not 
suffer from the resolution problem which plagues the FFT approach when only a 
short data record is available. It is recognized that the example presented is 
somewhat optimistic since neither noise nor the effects of some component 
error contributions were included in the simulation. The simulation does 
include non-ideal filter responses but does not incorporate the effects due to 
mixer imbalances. The plasma model accounts for reflections that occur at the 
vacuum-tile interface in addition to the reflections at the RCG boundaries. 
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Figure 3.5.1 Spectrum obtained from an autoregressive time series analysis of 
the group delay measurements made by the frequency-stepped 
DSBSC radar system. 

The intention of the example is to illuminate the fact that the technique of 
applying AR time series analysis to the frequency-stepped DSBSC 
measurements can be used to extract an accurate target distance estimate even 
when the SIR is significantly less than 0 dB. In our application of measuring 
electron plasma densities, a worst case Signal-to-Noise ratio (S/N) of 30 dB is 
expected at the inputs to the in-phase and quadrature detectors, and S/N ratios 
of 50-70 dB are expected for most measurements. Hence, it is not expected that 
noise levels will be large enough to significantly degrade performance in the 
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MRIS measurements. In other short-range radar problems, noise may be a 
significant problem and a more sophisticated spectral estimation approach such 
as an autoregressive moving average approach (ARMA) [Kay8t,88] may be 
required. 

It has been shown that autoregressive time series analysis can be used 
to measure accurately the period of the distance measurements made by the 
DSBSC system. The example presented shows that even when the SIR is 
significantly less than 0 dB, and the averaging technique fails as expected, the 
high resolution autoregressive estimator can provide an accurate distance 
estimate based on the simple two scatterer theory. 

3.6 Summary 

A new formulation of the frequency-stepped double sideband 
suppressed carrier radar system has been presented. With certain 
modifications, this chapter appeared in [Ybarra 2]. The newly derived system is 
an extension of the amplitude modulated (AM) CW technique [Nilssen] that was 
introduced initially in 1947 in [Riden]. The AM CW technique was abandoned 
because severely erroneous target range estimates were obtained when an 
interfering reflection was present. The main contribution of this chapter is the 
inclusion of the effect of an interfering reflection in the derivation of the DSBSC 
system. It is shown that the interference can actually be used as a reference 
allowing accurate target range estimates to be extracted from the radar data, 
independent of signal to interference ratio. 
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Simulation results were then presented in the context of measuring 
electron plasma densities that develop near the heat tiles of a space re-entry 
vehicle. The inhomogeneous media model derived in Chapter 2 was used to 
model EM propagation in the heat tile as well as the continuously tapered CFD 
predicted profile. The results illustrate that the DSBSC radar system is capable 
of accurate target range estimation even when the echo from the target is much 
weaker than the echo from the interferer. 
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Chapter 4 


The Frequency-Stepped Continuous Wave Radar System 

One of the most prevalent techniques for characterizing a network is to 
measure its response to several, different frequency sinusoidal signals. A few 
broad-class examples of the many applications which utilize this technique are 
measurement of transmission line discontinuities [Dworsky], bandwidth 
measurements [HP], and s-parameter determination of two-port networks 
[Collin]. There are two main reasons for the widespread use of this approach. 
First, frequency response measurements are easy to perform. All that is 
required is a source which generates sinusoidal signals (an oscillator) and a 
detector which measures the amplitude and phase of the response. Second, 
the frequency response of a linear, time-invariant network uniquely 
characterizes the system [Glisson]. Knowledge of a system's frequency 
response allows the computation of its response to any stimulus. These two 
factors combined have led to the production and use of literally thousands of 
network analyzers [HP], 

One of the primary objectives of a radar system is to determine the 
distance, or range, to one or several targets. Target range may be extracted 
from samples of the radar channel frequency response. This chapter presents 
an implementation of a radar system currently in use [Artech87] which employs 
a frequency-stepped CW approach to obtain measurements of the radar 
channel frequency response. There are two main objectives of this chapter. 
The first objective is to describe the mechanism by which frequency response 
samples are measured using the frequency-stepped CW approach. The 
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underlying principles of the stepped CW system are presented as the basis for 
the optimal signal processing algorithm derived and illustrated in Chapter 5. 
The second objective of this chapter is to illustrate the use of the 
inhomogeneous media model derived in Chapter 2 to simulate the 
measurement of non-uniform plasma densities using the frequency-stepped 
CW system. 

4.1. Introduction to the Frequency-Stepped CW System 

The physical data measured by the frequency-stepped CW radar system 
is a sequence of complex reflection coefficients as seen from the terminal plane 
to which the system is calibrated. The electric field reflection coefficient is 
defined by [Collin] 



c i (4.1.1) 

The reflection coefficient is, in general, a function of frequency and may be 
interpreted as samples of the channel frequency response H(jw). Frequency 
response for electric circuits is defined in general to be the ratio of the phasor 
output (voltage or current) to the phasor input (voltage or current). The 
analogous definition for EM waves is the ratio of the phasor output (electric field 
or magnetic field) to the phasor input (electric field or magnetic field). In the 

case of the reflection coefficient as frequency response, the output phasor 
quantity is the reflected electric field intensity E f and the input phasor quantity is 

the incident electric field E r 
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In the frequency-stepped CW system, a sinusoidal microwave or 

millimeterwave source (e.g. Gunn diode [Collins]) is used to generate the 
starting frequency f start of the measurement sequence. This carrier frequency is 

then mixed with another, lower frequency sinusoidal signal whose frequency is 
stepped in equally spaced increments Af. The frequency mixing performed in 

the transmitter is done in such a way as to produce an upper single sideband, 
frequency-stepped carrier whose frequency sequence is f start + iAf (i = 1 -n). 

The radar system is calibrated at a convenient terminal plane so that no 
reflection is measured at any measurement frequency when the system is 
terminated in a matched load. When a physical radar measurement is made, 
an interference pattern in the fields is created at the terminal plane to which the 
system is calibrated. This interference pattern is the resulting sum of the 
transmitted and reflected waves and is physically observed as the reflection 
coefficient at the terminal plane. The reflection coefficient is measured at each 
of the n discrete frequencies. Once the sweep is complete, a set of n 
measurements of the reflection coefficient, or frequency response, has been 
extracted. The remainder of this chapter examines this frequency response 
extraction process for a particular CW system implementation, and its 
performance is illustrated through simulations which employ the 
inhomogeneous media model derived in Chapter 2. The simulations utilize the 
standard IFFT approach to transforming the frequency domain data to the time 
domain, from which the target range is estimated. It is shown that the IFFT 
produces target range estimates which may contain a high degree of 
uncertainty. This range uncertainty produced by the standard IFFT approach is 
part of the motivation for the development of the new optimal signal processing 
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approach presented in Chapter 5. The main motivation for using the 
optimization approach is the poor resolution achievable using the IFFT. This is 
clearly demonstrated in Chapter 5 where experimental measurements made 
using a physical frequency-stepped CW radar system are analyzed using both 
the standard IFFT approach and the new optimal signal processing algorithm. 

4.2. The Frequency-Stepped CW Radar System Block Diagram 
and Principles of Operation 

The block diagram of the frequency-stepped CW system is shown in 
Figure 4.2.1, and is the configuration, except for the channel model, proposed 
by Electromagnetic Sciences [EMS]. The system operates as follows. Stairgen 
generates a staircase signal which drives the digitally controlled oscillator dco 
producing both an in-phase and quadrature sequence of sinusoids whose 
frequency is stepped from Af to 2Af to 3Af, through nAf. This frequency- 
stepped signal then modulates the in-phase and quadrature carrier. The mixed 
quadrature components are multiplied by -1 and summed with the mixed in- 
phase components producing an upper single sideband stepped carrier signal 
which is launched into the channel. The generation of the upper single 
sideband carrier is based on the trigonometric relation 

cos(27tf c t)cos(27rAf m t) - sin(27if c t)sin(27rAf m t) = cos[2jr(f c + Af m )t] (4.2. 1 ) 

The received signal is then beat against the carrier and filtered to retrieve the 
baseband frequency-stepped signal. 
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Figure 4.2.1 Block Diagram of the frequency-stepped CW radar system. 
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The coherent detector, which includes a lower single sideband generator, is 
then used to extract the in-phase and quadrature components of the frequency 
response data. Thus, the frequency-stepped CW system launches a sequence 
of n distinct unmodulated carriers into the channel, and the resulting induced 
gain and phase shift of the echo is measured by the receiver at each of the n 
frequencies. The system actually measures the in-phase and quadrature 
components of the echo which provide equivalent information. 

The amount of time spent on a single frequency step is called the step 
dwell time and its inverse is called the Step Repetition Frequency (SRF). The 
dwell time must be long enough so that the transients in the receiver filters have 
decayed sufficiently for steady-state measurement, and short enough so that the 
target does not move significantly during the n measurements. 

At the end of the n measurements, a set of n pairs of discrete frequency 
response data has been acquired. If the frequency change between steps is Af, 
then the total bandwidth over which the measurement is made is nAf. Using the 
standard approach, this finite bandwidth frequency response data is then 
converted into the time domain via an inverse Discrete Fast Fourier Transform 
(IFFT) to obtain an estimate of the channel impulse response. Through 
knowledge of the media dimensions and dielectric values through which the 
energy must propagate, the time domain impulse response may be converted 
into a distance domain impulse response. This distance domain impulse 
response usually possesses a peak value at the location of target. It will be 
shown in Chapter 5, however, that the IFFT is unreliable when the bandwidth is 
small and/or the target must be distinguished from an interfering reflection at 
approximately the same range as the target. 
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The next section presents the results of a computer simulation of the 
frequency-stepped CW system measuring the range to a specific electron 
plasma density using the cascaded slab model for the CFD plasma density 
profile. 

4.3 Simulation Results for the Frequency-Stepped CW System 

This section presents the results for a single simulation of the frequency- 
stepped CW system. It is intended to illustrate the use of the cascaded slab 
model for inhomogeneous media derived in Chapter 2, as well as some of the 
basic principles of using the IFFT to process the frequency domain data. The 
particular simulation presented employs the same channel and measurement 
bandwidth as was used to demonstrate the performance of the frequency- 
stepped DSBSC system derived in Chapter 3. Thus, a direct performance 
comparison of the two radar systems is possible by comparing the two target 
range measurements. Further results are presented for the stepped CW system 
in Chapter 5. 

Consider the CFD profile of figure 1 .2.1 scaled so that the critical electron 
density for 140 GHz is located at a standoff distance of 9 cm. In addition, let this 
plasma profile exist in the MRIS propagation path illustrated in figure 1.2.3. The 
time domain impulse response obtained using the IFFT on 64 measurements 
made by the frequency-stepped CW system over a bandwidth of 4 GHz, 
centered at 140 GHz is shown in figure 4.3.1. In this figure, there are actually 
two data sets plotted. The bar graph is the result of a 64 point IFFT of the 64 
complex frequency domain samples using a rectangular window. The 
continuous plot is the result of a 16, 384 point IFFT of the 64 frequency samples 
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padded with 32,704 zeroes (total IFFT length is 32,768 points). The zero- 
padding serves to interpolate between the time samples of the 64 point IFFT 
[Kay88] as illustrated in the figure. Using the dimensions of the vacuum, heat 
tile, and RCG coating, as well as their permittivities, the two-way travel time to 
the outside surface of the RCG coating is calculated to be 0.25 ns. Using c/2 as 
an approximation for the average group velocity in the plasma, a 9 cm standoff 
distance corresponds to a two-way travel time of 1.203 ns. Hence, the time 
delay for the reflection from the plasma turning point is 1 .453 ns. Examining the 
IFFT impulse response of figure 4.3.1 for this reflection it is evident that the exact 
reflection point is uncertain. The pulse resulting from the plasma has a peak 
value at a time delay of 1 .38 ns. Using this value as the time delay to the 
turning point, its range is calculated to be 8.475 cm standoff. The 
corresponding range error is only 5.25 mm. However, it is not clear that the 
peak in the impulse response should be at the turning point time delay because 
the pulse is not distinct. It is smeared in time due to the finite, 4 GHz bandwidth 
measurement window. The time domain impulse response is the result of the 
convolution of an impulse due to the reflection and a sine* pulse due to the 
rectangular window. The sine pulse is the IFFT of the rectangular window 
whose width is the measurement bandwidth of 4 GHz [Kay88] and is given 
analytically by 

F Gw) = = Ksinc(2jr*1 0 9 t) 

2tt* 1 0 t (4.3.1) 

where K is a constant. The time width of the main lobe of this sine function is 

sinx 

* sinc(x) is the sine-cardinal function and is defined here to be sinc(x) =■ — — . 
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(4.3.2) 


At ho = — 2- = 0.318 ns 

lobe 9 

27t*10 


This sine pulse is convolved with every impulse due to each reflection. If the 
simplifying assumption is made that there are only two reflections, one from the 
RCG and the other from the plasma, then the impulse response should consist 
of the sum of two sine functions each with a main lobe width of 0.318 ns. These 
sine pulses should be centered at 0.25 ns (RCG) and 1.453 ns (plasma). The 
sine pulses are clearly present in the impulse response of figure 4.3.1. 
However, these sine pulses interact. That is, the energy from the sidelobes of 
one sine function interferes with the other sine function. Consequently, the peak 
values of the two sine pulses are biased away from their values that would be 
obtained with no interference present. This is the analytical explanation for the 
resolution and accuracy deterioration in the ability of the IFFT to extract target 
range as the bandwidth of the measurement is decreased. As the bandwidth of 
the measurement is decreased, the sine functions spread out in time inversely 
proportional to the bandwidth. As the energy of the sine pulses spread out in 
time, the pulses interfere with one another, causing the locations of the peaks to 
change. Furthermore, if the bandwidth is reduced to a small enough value, the 
two sine pulses will merge into a single pulse and the distinction between two 
targets is lost. This phenomenon is carefully examined in the next chapter. 
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Figure 4.3.1 Impulse response obtained from 64 measurements of the MRIS 
plasma channel frequency response generated by the frequency- 
stepped CW radar system. The 64 point IFFT is shown along with 
a zero-padded 16,384 point IFFT. 


4.1. Summary 

This chapter has presented an implementation of a radar system 
currently in use [Artech87] which employs a frequency-stepped CW approach to 
obtain measurements of the radar channel frequency response. The 
implementation is a generalized system for producing CW measurements such 
as those measured by an HP 8510 network analyzer. The exact mechanism by 
which frequency response samples are measured using the frequency-stepped 
CW approach was described. The underlying principles of the stepped CW 


78 


system were presented as the basis for the optimal signal processing algorithm 
derived and illustrated in Chapter 5. Use is made of the inhomogeneous media 
model derived in Chapter 2 to simulate the measurement of non-uniform 
plasma densities using the frequency-stepped CW system. Target range 
uncertainty, inherent in IFFT processing due to the finite bandwidth of the 
measurement sweep, was clearly demonstrated in the results of the simulation. 
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Chapter 5 


Globally Optimal Signal Processing of the Frequency-Stepped CW 

Data 

The physical data measured by the stepped CW radar system is a 
sequence of complex reflection coefficients as seen from the aperture plane to 
which the system is calibrated. As was shown in Chapter 4, this sequence of 
measured reflection coefficients may be interpreted as samples of the channel 
frequency response H(jw). Given these measurements the objective is to 
determine the target range. The standard approach [HP] is to perform an IFFT 
which produces an estimate of the channel time domain impulse response. 
Peaks in the impulse response correspond to reflections and their time delay 
corresponds to the range to the reflection. The resolution of this approach is 
limited by the measurement bandwidth. Since the measurement bandwidth 
must be kept relatively narrow when the target is embedded in a frequency 
dependent inhomogeneous medium to avoid significant dispersion, the 
resolution offered by the IFFT approach may be unsatisfactory. 

In this chapter a new optimal signal processing algorithm is derived 
which maximizes the range resolution obtainable from any set of frequency- 
stepped CW measurements. The resolution limitations of the IFFT approach are 
demonstrated along with the enhancement offered by the new optimal 
processing algorithm using physical measurements made with an HP 8510 
network analyzer. Two different sets of measured data are analyzed. The first 
set is a sequence of reflection coefficients measured at the input port of the 
7mm Beatty Standard [Beatty]. The second set is a sequence of reflection 
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coefficients measured at the input port of a test apparatus for the MRIS 
instrument described in Chapter 1 . In order to illustrate the resolution 
enhancement achievable by optimal processing over the IFFT, the 
measurement bandwidth is reduced to the point that the IFFT approach clearly 
fails to resolve two distinct reflections that are known to exist. The optimal 
processing algorithm is then applied to the reduced bandwidth data and the 
results reveal the potential resolution achievable. 

5.1 The General Least Squares Optimization Problem 

Given a set of frequency response measurements generated by the 
frequency-stepped CW radar system, the objective is to determine the target 
range. This may be accomplished by assuming an underlying physical model 
for the channel, computing the frequency response of the model at the 
frequencies at which the measurements were taken, and finally minimizing the 
norm of the difference between the physical measurements and the values 
produced by the model. The analytical expression for the objective function to 
be minimized, denoted by J, is given by 

J = ||H(j(w 0 +iAw)) - H m G(w 0 +iAio))|| ( 5 .1 .1 ) 

where the H m (j(w 0 +iAio)) are the measured complex reflection coefficient pairs, 
the H(j(co 0 +iAw)) are the values of a theoretical model, || • || denotes vector norm 
[Golub], and i is the index number of each measurement and spans the integer 
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range from 1-n. J is a function of the parameters of the physical model yet to be 
specified. 

The physical propagation path which will be assumed contains a 
reflective target embedded in an inhomogeneous medium. The inhomogeneity 
may be either continuously tapered, as in the case of the non-uniform CFD 
plasma profile, or composed of distinct media layers each having a different 
permittivity, as in the case of the thermal protection heat tiles of the AFE 
spacecraft, or any combination of the two. The simplest realistic model for the 
system impulse response has the form 

h(t) = A 1 6(t-t 1 ) + A 2 5(t-t,) + A 3 5(t-t 3 ) + ... (5.1 .2) 

where the A ( are the reflection amplitudes and the tj are the time delays to the 
reflections. The corresponding frequency response of the system is obtained by 
taking the Fourier Transform of (5.1 .2) which produces 

H(j(o) = + Age' 1 ^ + A 3 e' jl0, =» + ... (51 3 ) 

Samples of this frequency frequency response form the model values 
comprising H(j((o 0 +iAw)) in (5.1.1). 

H(j(w 0 +iAto)) = + A 2 e‘ K “o +iA “> t 2 + A 3 e' j(oJ o +iAo3 > t 3 + ... ^ 

Therefore, using the assumed model for the system impulse response (5.1 .2), 
the objective function J (5.1.1) is a function of the reflection amplitudes and 
delays. 

The type of vector norm implemented in (5.1.1) must be judisciously 
chosen. A proper choice of norm will incorporate all the measured data and 
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present a tractable mathematical problem for extracting the optimal values for 
the amplitudes and delays of the model. A reasonable choice is the 2-norm 
since the objective function J will be quadratic in the amplitude parameters and 
a standard least squares approach [Golub, Lawson] can be used to extract 
those parameters. Using the 2-norm as the metric in (5.1.1), the general least 
squares optimization problem becomes finding the set of amplitudes and delays 
in the model which minimize 

J = £ |H(j(to 0 +iA(o)) - H m (j(w 0 +iAw))| 2 

hi (5.1.5) 

The summation index / in (5.1 .5) ranges over the n measured data pairs. In the 
next section a derivation is presented for determining the optimal reflection 
amplitudes and delays. 

5.2 Derivation of the Optimal Reflection Amplitudes 

The objective is to find the global minimum of the performance metric 
(5.1 .5) with respect to the amplitude-delay pairs A jf tj in the model (5.1 .2-3). This 

objective is achieved by decomposing the performance metric (5.1 .5) into a sum 
of linear and non-linear components. The linear portion of the objective is a set 
of simultaneous linear equations whose solution yields the optimal amplitudes. 
The metric (5.1 .5) is quadratic in the amplitude parameters Aj. This allows the Aj 

to be calculated in closed form by solving the least squares problem 

Fa * f (5.2.1) 
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The complex (n x N) matrix F is composed of frequency response estimates 
based on the model (5.1 .4) and is given by 


cos(u3 1 1^) - j sin(ti3 1 t N ) 
COS^tfg) - 1 Sin(oo 2 t|\i) 

cos(oo n t N ) - j sin(u n t N ) 

(5.2.2) 

where n is the number of frequency measurements and N is the number of A jf tj 
pairs to be extracted. The complex vector f is composed of the physical 
frequency response measurements and is given by 

f = IH M (w 1 ) H m (w 2 ) Hjcog) ••• H M (w n )] T (5.2.3) 

The vector a in (5.2.1) contains the set of amplitudes Aj to be determined 

a =[A 1 Aj, A 3 - A,/ ( 5 2.4) 

Solution of the least squares problem requires assuming a set of delays tj. It is 
shown later that by "scanning" a sequence of amplitude optimized objective 
function values for certain delay combinations, it is possible to find the global 
minimum of (5.1.5). Finding the global, or absolute minimum of (5.1.5) 
determines the optimal set of reflection amplitude-delay pairs. If there exists 
any a priori knowledge about the target range such as upper and/or lower 
bounds, the range of the time delays over which the objective function is 
scanned can be limited and the amount of computation necessary to find the 
global minimum reduced. 


F = 


cos(w 1 t 1 )-jsin(o5 1 t 1 ) cos(u).,t 2 ) - jsin(w.,t 2 ) 
cos(w 2 t 1 ) - i sin(oo 2 t 1 ) cos(u 2 t 2 ) - j sin(to 2 t 2 ) 

cosfapt.,) - j Sin(i0 n t 1 ) cos(u5 n t 2 ) - j sin(w n t 2 ) 
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The solution to (5.2.1) will yield a set of complex amplitudes because the 
model (5. 1.2-4) will not fit the physical data exactly. However, in the derivation 
which follows, the amplitude vector a is constrained to be real. The same 
approach that is used to solve the standard real LS problem using the normal 
equation approach [Golub, Lawson] can be extended to solve this problem. 
The solution to (5.2.1) is found by minimizing the real scalar that results from the 


following squared inner product. 


y = II (F, + j F^a - (f r + j t,)ll| 

(5.2.5) 

y = (F,a - f r ) T (F,a - «,) + (F,a - f,) T (F ( a - »,) 

(5.2.6) 

y = a T [F>, + F T F,]a - a T [F*f, + F T f,] - [f>, + ffaa +!>, + »7 f , 

(5.2.7) 

In order to minimize y, the gradient of (5.2.7) is taken with respect to the 
amplitude vector a and set equal to the zero vector. 

|| =0 = 2[FjF, + F> l la-2[F\ + F[f i ] 

(5.2.8) 

The result of the minimization of (5.2.5) is a real, square set of linear equations 
which can be solved by using any linear equation solving technique. 

IF^F, + F^F.]a = F^‘ r + F^V 

(5.2.9) 


However, in certain cases, the symmetric matrix premultiplying a is extremely 
ill-conditioned [Golub, Lawson]. In order to see when such cases arise let 
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(5.2.10) 


H=F T F , + F> 


i 1 i 


and examine the case of a two reflection problem for simplicity. 

£ cos 1 2 3 (Wjt ? ) + sin 2 (c0jt 2 ) = n £ cos(«jt.|)cos(i0jt 2 ) +sin(Wjt 1 )sin(w j t 2 ) 
i =1 1=1 

l cosfWjt^cosfWjtg) + sintWjt^sintWjtg) .X cos 2 (oi s t 2 ) + sin 2 (to s t 2 ) =n 

( 5 . 2 . 11 ) 


H 


2x2 


One obvious case that will cause H to be ill-conditioned occurs when the values 
of t.j and t 2 are nearly equal. When the values of t 1 and t 2 are exactly equal, H 

becomes singular. In order to deal with this problem, the singular value 
decomposition technique [Golub, Press] is used to solve (5.2.9). When any 
singular value is found to be less than 0.01 , its value is set equal to zero before 
the final computation of a is performed. This extracts the minimum norm 
solution to (5.2.9), which is precisely what is sought. 

The original LS optimization problem (5.1.5) has been decomposed into 
a linear problem whose solution yields the amplitudes Aj, and a highly non 
linear problem whose solution, it will be shown, yields the delays tj. It is 
important to recognize that the solution of (5.2.9) for the amplitudes Aj requires 

that a set of delays tj be utilized in forming the entries of the matrix H. 

The following optimization procedure may be followed for extracting the 
Aj, tj pairs: 

1) Choose the number N of Aj, tj pairs to be determined. 

2) Assume a sequence of tj values and solve (5.2.9) for the optimal Aj . 

3) Search for the global minimum of (5.1 .5) using the amplitudes found from 
the solution of (5.2.9). 
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It will be shown in Section 5.4 that the geometric interpretation of the LS 
objective function J, subject to the solution of (5.2.9), is a multiple minima 
surface with grooves aligned with the delay axes. The grooves in this surface 
can be searched sequentially and the global minimum found with little 
computation. The global minimum corresponds to the optimum reflection 
amplitudes and delays and hence to target range. The next section presents 
frequency-stepped CW measurements made using an HP 8510 network 
analyzer which will then be analyzed using the global optimization technique in 
section 5.4. 

5.3 Experimental Measurements 


This section presents frequency-stepped CW measurements made using 
an HP 8510 network analyzer in two different experimental set-ups. First, 
measurements of the reflections in an air-line coaxial cable section, the Beatty 
Standard, are presented. Then, measured data from an MRIS test apparatus 
including a simulated heat tile is examined. The process of extracting the 
reflection time delays using the standard IFFT approach, as used by the HP 
8510 network analyzer internal computer, is illustrated for both sets of 
experimental data. 

Consider the 7 mm Beatty Standard [Beatty, HP] of the HP-8510 network 
analyzer terminated in a matched load and its theoretical bounce diagram as 
shown in figure 5.3.1 . The Beatty Standard presents two step-discontinuities in 
the characteristic impedance of the line. These discontinuities are produced by 


an abrupt change in the outer diameter of the inner conductor. The bounce 
diagram shows the theoretical amplitude and delay of each reflection for 

measurement calibrated to the input port of the Beatty Standard. 
Measurements of the reflection coefficient s n were taken over two different 

bandwidths using an HP8510 network analyzer. 

The first set of measurements were taken over the band 45 MHz - 18 GHz 
using 801 evenly spaced frequencies. The resulting IFFT is shown in figure 
5.3.2. A rectangular window [Harris] was utilized to maximize the resolution 
between the reflections. The 18 GHz bandwidth is large enough to clearly 
distinguish the reflections. In addition, the amplitude and delay of the first three 
reflections is very accurate. The second set of measurements were taken over 
the band 4 GHz - 6 GHz using 801 evenly spaced frequencies, and the resulting 
IFFT is shown in figure 5.3.3. 
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Figure 5.3.1 7 mm coaxial Beatty Standard terminated in a matched load and 

the resulting bounce diagram for a unit impulse stimulus. 






|h(t)l 


Figure 5.3.2 


Figure 5.3.3 



Inverse FFT of 801 measurements of the reflection coefficient at 
the input port of the 7 mm Beatty Standard shown in figure 5.3.1 . 
The bandwidth used in this measurement is 18 GHz. 


|h(t)| 



Inverse FFT of 801 measurements of the reflection coefficient at 
the input port of the 7 mm Beatty Standard shown in figure 5.3.1 . 
The bandwidth used in this measurement is 2 GHz. 
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The 2 GHz bandwidth is so small that the two reflections (83 ps, 583 ps) have 
merged into a single pulse. In this case the IFFT cannot resolve the two 
reflections. In addition, the location of the peak in the impulse response is 
biased [Kay88] away from the correct location of either of the two reflections. 
Hence, when two reflections are so close in time that the IFFT cannot resolve 
them, accuracy is lost as well. 

In order to test the applicability of the optimization approach to detecting 
the turning point range for the MRIS instrument, a test fixture was constructed at 
NASA LaRC. This fixture is shown in figure 5.3.4 along with its physical 
connection to the network analyzer via a rectangular to circular waveguide 
transition. The network analyzer is calibrated to the input port of the waveguide 
as shown in the diagram. The simulated heat tile is composed of a 0.25 thick 
layer of silicon dioxide, or quartz glass (e r = 3.78), and a 0.5" thick layer of 
polyurethane foam (e r = 1 .4). A 0.5" thick aluminum plate is used for the 

aperture ground plane. The target is formed by the presence of another 0.5" 
thick aluminum plate parallel to the aperture plane of the circular horn antenna. 
Target range is set by moving mobile target. 

In the first set of measurements, the target was set at a distance of 12" 
from the outer surface of the quartz glass. The two primary reflections will be 
due to the high dielectric constant quartz glass and the metal plate target. 801 
measurements of the reflection coefficient were made from 14-18 GHz, 
spanning a total bandwidth of 4 GHz. The IFFT of the measured radar data is 
shown in figure 5.3.5. The first three reflections are clearly resolved. The time 
delay corresponding to the first reflection, due to the quartz glass interferer, is 
difficult to calculate accurately because the velocity of propagation throughout 
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the waveguide and antenna is not known. The velocity of energy propagation, 
or group velocity, within a waveguide is frequency dependent and is given by 
[Collin] 


v 


g 



(5.3.1) 


where c is the free space velocity of light, f c is the cut-off frequency of the 
waveguide, and f is the frequency of propagation. The rectangular waveguide 
in figure 5.3.4 has the dimensions 1 .58 cm by 7.9 mm. This geometry produces 
a cut-off frequency of 9.49 GHz [Collin]. In order to determine the average 
group velocity, Vg must be integrated over the measurement bandwidth. 


v = — £ 


aV f 2" f 1 


df 


'1 


(5.3.2) 


The antiderivative is of standard form and from [CRC] is given by 


f 2 " f 1 


1-1^1 df = r - £ 


f 2 * f 1 




f sec' 1 !— 
c c |f 

c 


(5.3.3) 


Using the frequency span of 14-18 GHz in (5.3.2) and a cut-off frequency of 9.49 
GHz, the average velocity of propagation is found to be 0.801c within the 
rectangular waveguide. If the entire length of waveguide in figure 5.3.2 were 
rectangular, the two-way travel time to the center of the quartz glass would be 
3.71 ns. The time delay to the first major reflection is approximately 3.45 ns as 
determined from the IFFT shown in figure 5.3.5. The time delay calculation 
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above, which assumes that the entire waveguide and antenna system is 
composed of a uniform rectangular waveguide, is intended only as an example 
of an approach which may be used to determine the average propagation delay 
through the waveguide system. The calculation of the average group velocity 
within the waveguide structure is actually more complicated due to the 
rectangular to circular transition as well as the flare in the circular horn antenna. 
The rectangular to circular waveguide transition is designed to suppress the 

creation of multiple modes. A gradual transition from rectangular to circular 
geometry allows the energy in the dominant mode (TE 10 ) of the rectangular 

guide to propagate as energy in the dominant mode (TE^) of the circular guide 
[Collin], However, the circular horn stimulates the creation of several TE and 
TM modes making the calculation of average group velocity more complicated 
than assuming single mode transmission. By carefully examining the geometry 
of the waveguide and circular horn antenna, the average propagation velocity 
may be calculated leading to a very accurate prediction of the time delay to the 
quartz glass. 
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Figure 5.3.5 Normalized Inverse FFT of 801 measurements of the reflection 
coefficient at the input port of the MRIS test apparatus shown in 
figure 5.3.4. The bandwidth used in this measurement is 4 GHz, 
and the target range is 12". 

However, the prediction of the time delay to the interfering reflection from the 
quartz glass is not the primary interest here. The primary concern is target 
range resolution and accuracy. Detecting accurately the time delay from the 
reflection due to the quartz glass to the metal plate target is the objective. 

The IFFT shown in figure 5.3.5 indicates there is a delay of 1.97 ns from 
the quartz glass to the metal plate. Using the free space velocity of propagation 
for the two-way distance of 24" in air requires 2.032 ns. Since each reflection 
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from from the quartz glass is assumed to be a single impulse, two unresolved 
impulses resulting from the front and back surfaces, an additional delay of 41 ps 
is included to account for the 0.25" of travel through the glass. Thus, the target 
range is defined to be the distance from the center of the glass to the metal 
plate. The total theoretical two-way travel time is 2.073 ns. A time delay 
estimate of 1 .97 ns is produced by the IFFT. Thus, the error in the time delay 
estimate produced by the IFFT using 4 GHz of bandwidth is 103 ps and 
corresponds to a range estimate error of approximately 0.608" (1 .545 cm). 

In the second set of measurements, the target was set at a distance of 
2.25" from the outer surface of the quartz glass. 100 measurements of the 
reflection coefficient were made from 14-14.5 GHz, spanning a total bandwidth 
of only 500 MHz. The IFFT of the measured radar data is shown in figure 5.3.6. 
The narrow bandwidth spanned by the measurement sequence is so small that 
several reflections have merged into a single lobe of the IFFT. The 2.25 
distance from the outer surface of the quartz glass to the metal plate target 
corresponds to 0.381 ns of two-way travel time. Including 41 ps of delay within 
the quartz glass results in a total theoretical two-way travel time of 0.422 ns. 
Since the metal plate target is so close to the tile, there are many reflections 
between the quartz glass and the aluminum plate. These reflections are 
separated by 0.422 ns in time. The main lobe in the IFFT spans approximately 
2 ns, enough time to "swallow" nearly 3 of these multibounce reflections. The 
time delay corresponding to the second major peak in the IFFT of figure 5.3.6 is 
4.8 ns and has no value in assessing the range to the target. This is the best 
performance the IFFT has to offer for determining the distance to the target 
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under the bandwidth constraint of 500 MHz with only 100 frequency 
measurements. 


|h(t)| 



Figure 5.3.6 Normalized Inverse FFT of 100 measurements of the reflection 
coefficient at the input port of the MRIS test apparatus shown in 
figure 5.3.4. The bandwidth used in this measurement is 500 
MHz, and the target range is 2.25". 


5.4 Optimal Signal Processing of the Experimental Measurements 


This chapter presents the results of applying the globally optimal signal 
processing technique derived in section 5.2 on physical measurements. The 
processing technique is globally optimal in the sense that the minimized 
objective function (5.1.5) has many local minima and a search procedure is 
used to find the absolute minimum. This global minimum determines the best 
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set of reflection amplitudes and delays given the frequency data and assuming 
that the radar channel has distinctly reflective scatterers. 

Consider the objective function (5.1 .5) for two assumed delays, one 
scanned from 0-0.8 ns along one delay dimension and the other from 0.06-0.1 
ns using the Beatty Standard data and 18 GHz bandwidth. This function may 

be interpreted as a surface above the "delay" plane as shown in figure 5.4.1 . It 
is very important to recognize that every value of J(tj, tj) is the result of solving 
(5.2.9) using each delay pair t jf tj in the domain shown. The surface shown in 

figure 5.4.1 has a single global minimum which corresponds the optimum delay 
pair (t 1 = 83 ps, t 2 = 583 ps). The grooves in the objective function extend along 

the delay axes. The deeper groove at 0.083 ns is due to the larger amplitude 
reflection caused by the first discontinuity of the Beatty Standard coaxial line 
(amplitude = -0.3333). The groove at 0.583 ns is due to the slightly smaller 
amplitude reflection from the second discontinuity (amplitude = 0.297). The 
surface has multiple local minima, hence a simple descent algorithm will fail, in 
general, to locate the global minimum. One approach that will be demonstrated 
here is to use a grid search algorithm on the objective function. Grid search is a 
method of finding the global minimum by finding the least value of J along one 
delay dimension while holding the other delay fixed. Then, once the least value 
of J is found in that dimension, the corresponding delay is held fixed, and the 
search is performed along the other dimension. By performing multiple 
iterations in this fashion, the global minimum is found once the delays are no 

longer updated. The global minimum for the 18 GHz case (figure 5.4.1) was 
found in two iterations to be t 1 = 0.082 ns, t 2 = 0.588 ns, very close to both the 

theoretical delays as well as the estimates produced by the IFFT approach. 
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The more interesting case occurs when the bandwidth is reduced to 2 GHz and 
the IFFT is useless for determining the range to either of the first two reflections. 
Applying the optimal signal processing algorithm to the frequency data for this 
case results in the objective function illustrated in figure 5.4.2. There still exist 
grooves in the objective function. However there are several grooves whose 
minima lie between the values 0.083 ns and 0.583 ns. There does exist a 
single global minimum in this domain and the delay values obtained using grid 
search to find it are .082 ns and 0.587 ns, which is more accurate than the delay 
estimates produced by the IFFT for 18 GHz bandwidth. The impulse response 
generated by the optimization approach is shown in figure 5.4.3. This result 
alone shows the resolution enhancement offered by the optimization approach 
over the standard IFFT as used by the HP 8510 network analyzer. 
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Figure 5.4.2 Amplitude optimized objective function constructed from 801 
reflection coefficient measurements of the Beatty Standard 
(figure 5.3.1) using a bandwidth of 2 GHz. 
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Figure 5.4.3 Impulse response of the Beatty Standard obtained using grid 
search (3 iterations) on the amplitude-optimized objective 
function resulting from 2 GHz bandwidth HP8510 
measurements. 


A more impressive set of results is obtained from the radar 
measurements of the MRIS test apparatus of figure 5.3.4. The objective 
function surface for 801 measurements over a 4 GHz bandwidth is shown in 
figure 5.4.4. This surface clearly shows grooves along the delay axes. The 
domain of this plot was chosen to display the symmetry of the objective surface. 
As shown, there are two global minima one with t 1 = 3.47 ns, t 2 = 5.43 ns and 
the other with t 1 = 5.43 ns, t 2 = 3.47 ns, which are exactly the same solution to 

the problem. Now consider the case of only 100 measurements spanning a 
bandwidth of only 500 MHz. The IFFT produced the meaningless result shown 
in figure 5.3.6. The optimization approach produces the objective surface 
shown in figure 5.4.5. 
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Fiaure 545 Amplitude optimized objective function constructed from 801 
9 reflection coefficient measurements of the MRIS test apparatus 

(figure 5.3.4) using a bandwidth of 500 MHz. Target is at a 
distance of 2.25" from the outer surface of the quartz glass. 


This objective function surface may appear to lack structure. However, the 
surface has grooves running along the delay axes which can be scanned 
successfully using the grid search algorithm. Evidence of these grooves may 
be seen by carefully examining the right half of the surface. The ridges In the 
surface indicate the presence of the grooves. The surface also possesses a 

macroscopic bowl shape, albeit very "noisy," in the domain shown. There are 
two global minima, either of which produces the delay estimates t 1 = 3.53 ns, t 2 

= 3.95 ns. These delay values were found by polishing the result of a first-pass 
grid search with an additional grid search. In the first pass, the domain 

scanned was 1-6 ns in each delay dimension using 50 ps increments. After 16 
iterations of grid search, the resulting time delays are t 1 = 3.55 ns, t 2 = 3.95 ns. 

Then, a second pass was performed over the domain 3.5-3.6 ns in the first 
delay dimension and 3.9-4. 0 ns in the second delay dimension using 1 ps 
increments. After 12 iterations of grid search in the second pass, the resulting 
time delays are t 1 = 3.53 ns, t 2 = 3.95 ns. The impulse response produced by 

the optimization technique is shown in figure 5.4.6. 

One interesting observation that should be made is the that the estimate 
of the time delay to the quartz glass is larger (3.53 ns) for the 500 MHz case 
than for the 4 GHz case (3.45 ns, see figure 5.3.5). Thus, it appears that the 
glass has moved in range. The cause of the 80 ps additional delay is due to the 
fact that the frequency scan in the 500 MHz case is from 14 GHz to 14.5 GHz, a 
band whose average frequency is less than the average frequency for the 4 
GHz case which scanned from 14-18 GHz. Hence, from (5.3.2) the average 
group velocity is less for the 500 MHz sweep than for the 4 GHz sweep, and a 
longer time delay to the glass is expected. 
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The most important measurement is the difference in time delay between 
the two reflections. This delay difference is found to be t 2 -t 1 = 0.42 ns. The 

theoretical value was calculated to be 0.422 ns. Hence, a delay error of only 
2 ps is produced by the optimization approach in this case. This 2 ps error in 
the delay estimate corresponds to approximately 0.3 mm. This result clearly 
shows the resolution capability of the optimization approach. The IFFT for this 
measurement case (figure 5.3.6) was totally useless for determining the range 
to any reflection while the optimization approach has distinctly and accurately 
resolved both the reflection from the metal plate target as well as the reflection 
from the quartz glass. 


|h(t)| 



time (ns) 


Figure 5.4.6 Normalized Impulse response of the MRIS test apparatus of 
figure 5.3.4 obtained using grid search (16 iterations, 14 
iterations polished) on the amplitude-optimized objective 
function resulting from 500 MHz bandwidth HP8510 
measurements. Target range is 2.25". 
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The results for the IFFT and optimization approaches for both the Beatty 
Standard and MRIS experiments are presented in the tables below. 


Table 5.4.1 Comparison of reflection delay estimates produced by the IFFT 
and optimization approaches for the Beatty Standard. 


BW =18 GHz, n=801 measurements 

Theoretical 

IEEE 

Optimization 

Time delay to the first reflection: 

83 ps 

82 ps 

82 ps 

Time delay to the second reflection: 

583 ps ' 

585 ps 

588 ps 

BW = 2 GHz, n=801 measurements 

Theoretical 

IEEE 

Optimization 

Time delay to the first reflection: 

83 ps 

N/A* 

82 ps 

Time delay to the second reflection: 

583 ps 

N/A* 

587 ps 

Table 5.4.2 Comparison of target delay estimates produced by the IFFT and 
optimization approaches for the MRIS test apparatus. The 
distance from the outer surface of the quartz glass to the metal 
plate target is 12" for the 4 GHz case and 2.25 for the 500 MHz 
case. All time delays shown are two-way travel times from the 
center of the quartz glass to the metal plate. 

BW = 4 GHz, n=801 measurements 

Theoretical 

1EEI 

Optimization 

Time delay from the quartz glass 
to the metal plate target: 

2.073 ns 

1.97 ns 

1.96 ns 

BW = 500 MHz, n=100 measurements 

Theoretical 

IEEE 

Optimization 

Time delay from the quartz glass 
to the metal plate target: 

0.422 ns 

N/A* 

0.420 ns 


* It is not clear from the IFFT results that there are two reflections. Even if 

assumed, the time delays to the peaks in the IFFT results are erroneous to the point of beng 

meaningless. 
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An existing parametric technique called the Matrix Pencil Method [Hua, 
Maricevic] was designed to recover amplitude-delay pairs from frequency 
domain data. The Matrix Pencil Method has been shown to perform 
significantly better than the IFFT approach. The results presented in this 
chapter for the optimization approach applied to the Beatty Standard are 
comparable to those presented for the Matrix Pencil Method [Maricevic]. 
However, the optimization approach is much more general. For example, time- 
gating [Artech84,87] can be implemented easily in the optimization approach. 
This is equivalent to searching only over certain regions in the delay plane. No 
physical constraints such as time-gating can be implemented in the Matrix 
Pencil Method without a significant modification in its derivation. 

5.5 Summary 

This chapter has presented the derivation and demonstration of a new 
algorithm for processing radar data produced by any frequency-stepped CW 
system. The derivation is based on minimizing the two-norm of the difference 
between the sequence of measured reflection coefficients and those produced 
by a model which assumes that the echoes from the radar channel are impulses 
in the time domain. A performance comparison between the IFFT and optimal 
processing techniques was presented in the context of two physical 
measurement sets. In the first set, the Beatty standard was terminated in a 
matched load and driven by an HP 8510 network analyzer. The measurement 
bandwidth was reduced from 18 GHz to 2 GHz and it was shown that while the 
IFFT produced meaningless reflection range estimates in the 2 GHz case, the 
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optimal processing approach was able to clearly resolve the reflections and the 
resulting reflection range estimates were more accurate using 2 GHz bandwidth 
than those produced by the IFFT using 18 GHz. In the second set of 
measurements, an MRIS test apparatus was also driven by an HP 8510 network 
analyzer. In this set of measurements the bandwidth, target range and the 
number of measurements were all reduced. First, the metal plate target was set 
at a range of 12" and 801 measurements spanning a total bandwidth of 4 GHz 
were analyzed using both the IFFT and optimal processing algorithms. The 
target range estimates resulting from both methods were of equal accuracy. 
Then the target was set to a range of 2.25", the number of measurements was 
reduced to 100 and the bandwidth was reduced to only 500 MHz. The range 
resolution and accuracy deterioration produced by IFFT processing was clearly 
evident. In this case, the IFFT could not distinguish the presence of two 
reflections nor could it determine accurately the range to either the interfering 
reflection or the metal plate target. However, the optimal processing approach 
could not only resolve the two reflections but could estimate both of their ranges 
with high accuracy. This result demonstrates the resolution enhancement 
offered by the optimal processing approach over IFFT processing. 
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Chapter 6 

Conclusions and Future Research 


This dissertation provides three original contributions to the field of 
microwave-millimeterwave radar. A model for electromagnetic planewave 
propagation in inhomogeneous media was derived in chapter 2 and its 
performance illustrated in the context of modelling EM propagation in non- 
uniform electron plasmas. This model is based on the theory of transmission 
lines and represents an inhomogeneous medium by thin cascaded slabs, each 
of which is homogeneous but whose combined structure is inhomogeneous. It 
has been shown that as the thickness of each slab is reduced, the response 
from the model asymptotically approaches the true physical response. The 
model can be implemented very efficiently using advanced computer 
programming techniques such as recursive programming which allow the use 
of hundreds of thousands of slabs while requiring little computation time. 

A new formulation for the frequency-stepped double sideband 
suppressed carrier radar system was derived in Chapter 3. This reformulation 
included the effect of an interfering reflection. Simulation results were 
presented which illustrate that the DSBSC radar system is capable of accurate 
target range estimation even when the echo from the target is much weaker 
than the echo from the interferer. 

A new globally optimal signal processing algorithm was derived in 
Chapter 5 for processing radar data produced by the frequency-stepped CW 
system. A non-linear least squares objective function is formulated, and its 
absolute minimum yields both the delay times and amplitudes of the reflections 


110 


from the target and any reflective interferer. The performance of the algorithm is 
demonstrated using physical radar data and it is shown to produce accurate 
and highly resolved target range estimates even when the bandwidth is 
reduced well beyond the point where the standard IFFT approach fails. 

Several simulation results have been presented in the context of 
measuring non-uniform electron plasma densities that develop near the heat 
tiles of a space re-entry vehicle. Unfortunately, an appropriate plasma profile 
cannot be produced for ground testing. Hence, with the exception of performing 
an actual atmospheric aeropass, their measurement can only be performed by 
computer simulation. For this reason, the inhomogeneous media model 
derived in this dissertation, which can simulate accurately the propagation of 
planewaves in an electron plasma, is especially significant. 

Several research problems are generated by the work presented in this 
dissertation. The inhomogeneous media model should be extended to include 
oblique incidence which will require the effects of polarization to be included as 
well. In [Pozar] the transmission line analogy for EM planewave propagation 
through dielectric slabs is extended to include oblique incidence. Hence, the 
framework for extending the cascaded slab model to include oblique incidence 
already exists. Another extension of interest would be the inclusion of antenna 
radiation pattern [Jordan, Silver]. The MRIS system described in Chapter 1 
utilizes a circular horn antenna mounted behind the heat tiles of the AFE 
spacecraft. The radiation transmitted into the propagation path is approximated 
by planewaves incident normally on each of the media layers. However, 
because the antenna is so close to these media layers, the radiation is actually 
composed of spherical waves and therefore much of the radiation is obliquely 
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incident on the layer boundaries. The intensity of radiation as a function of 
angle constitutes the antenna radiation pattern and must be included in the 
inhomogeneous media model for a truly accurate representation of the physical 
measurement. 

In the derivation of the DSBSC system in Chapter 2, as well as the 
stepped CW system in Chapter 3, it was tacitly assumed that both the target and 
interferer are stationary. If either the target or interferer move significantly 
during the time in which a measurement is made, the received echoes will 
contain a Doppler shift. The effect of target motion on the measurements made 
by these radar systems remains to be explored. Standard radar techniques 
such as pulsed radar, extract target velocity from the Doppler shifted echo. It is 
expected that target velocity estimation may be achieved by both frequency- 
stepped radars examined in this dissertation. However, this remains to be 
shown. 

The DSBSC system was derived for exactly two reflections, one from the 
target and one from the interferer. This derivation provides a significant 
improvement in the applicability of the DSB technique. However, the technique 
should be generalized to include the effects of several reflections. 

The derivation of the optimal signal processing algorithm in Chapter 5 for 
the stepped CW system made no assumption about the number of reflections. 
However, the simulations presented examined only the case of extracting two 
reflections from the radar data. For the case of two reflections, it was 
demonstrated that a simple grid search procedure can be used to find the 
global minimum of the objective function. The optimization process should be 
investigated for an arbitrary number of reflections, it is likely that the simple grid 
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search technique will fail in the more general case. A more sophisticated 
search technique, such as simulated annealing [Bev], may be required. 

There may be applications for the two radar systems presented in this 
dissertation for which the propagation path induces significant noise in the 
radar measurements. This issue has not been addressed and remains as 
future research. 
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Chapter 8 

Appendix: Derivation of Equation 3.2.14 


r . 2(a-M /Vc2-d 2 tan(|) 

( a+beos(x) dy _ fax. + 3 — stl tan ' 1 

' ' J c+dcos(x) d f c 2 . d 2 


c + d 


(8.1) 


Proof: 


A substitution of variables will be made which will transform the 
trigonometric polynomial integrand of F(x) into a rational polynomial. Let 

u = tan (f) (8.2a) 


This produces 


x = 2 tan’ 1 u dx = — - - du 

1 + u 2 


(8.2b) 


and 


cos(x) = cos(2tan‘ 1 u) = 1 - 2sin 2 (2tan' 1 u) = 1 - 2sin 2 (sin' 1 (^=M=r| (g3) 


cos(x) = 1 


_2uL = 

1+u 2 1+u 2 


(8-4) 


F(x) may be decomposed into two simpler functions F^x) and F 2 (x). 


a+bcos(x) dx _ a | dx 
c+dcos(x) I c+dcos(x) 



c ° s(x) , dx = F,(x) + F 2 (x) 
c+dcos(x) 


(8.5) 
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Transforming the function F^x) into the function F^u) using the transformation 
(8.2) yields 



Recognizing the integrand of (8.7) as a having the familiar derivative relation 

-d-tan'W = — du 

du VP/ p +u 2 (8.8) 


allows F.,(u) to be expressed as 


Simply replacing 


Fv(u) = ^ A / r ^4tan' 1 (u A /^J] 
n ; c-dVc+d \ Vc+di 


the variable u with tan(x/2) gives 



(8.9) 


( 8 . 10 ) 
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Now transforming the function F 2 (x) into the function F 2 (u) using the 
transformation (8.2) yields 


F 2 (u) = b 


J 


r^j2 

_1±lj 2 2__du 

H ^) 1+u2 


= 2b f 1=U? 

J (c(1 +u3)+d(1 -u 2 )X 1 +u3) 


du 


( 8 . 11 ) 


Fs(u)= ^ 

c-d 


1-U 2 


|c+d + U 2 j(l+U 2 ) 


du 


( 8 . 12 ) 


The integrand of (8.12) may be decomposed into two simpler functions using 
partial fraction expansion giving 


f 


p2(U) = I 




B 

(l + u 2 ) 


du 


(8.13) 


where the constants A and B are easily determined to be 



Inserting these coefficients into (8.13) gives 


r 



(8.14) 
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The first anti-derivative in (8.14) is identical to (8.7) and the second is easily 
evaluated using the derivative relation (8.8) yielding 


F 2 (u) = 



+ 


^tan'V) 

d 


(8.15) 


Once again substituting u=tan(x/2) 






/ \c 2 -d 2 tan(£)\ 


(8.16) 


Combining F^x) with F 2 (x) completes the proof. 
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